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ABSTRACT 

A  nonlinear  representation  for  the  rapid-pressure  correlation  appearing  in  the  Reynolds  stress  equa¬ 
tions,  consistent  with  the  Taylor- Proudman  theorem,  is  presented.  The  representation  insures  that 
the  modeled  second-order  equations  are  frame-invariant  with  respect  to  rotation  when  the  flow  is 
two-dimensional  in  planes  perpendicular  to  the  axis  of  rotation.  The  representation  satisfies  real¬ 
izability  in  a  new  way:  a  special  ansatz  is  used  to  obtain,  analytically,  the  values  of  coefficients 
valid  away  from  the  realizability  limit:  the  model  coefficients  are  functions  of  the  state  of  the  tur¬ 
bulence  that  are  valid  for  all  states  of  the  mechanical  turbulence  attaining  their  constant  limiting 
values  only  when  the  limit  state  is  achieved.  Utilization  of  all  the  mathematical  constraints  are 
not  enough  to  specify  all  the  coefficients  in  the  model.  The  unspecified  coefficients  appear  as  free 
parameters  which  are  used  to  insure  that  the  representation  is  asymptotically  consistent  with  the 
known  equilibrium  states  of  a  homogeneous  sheared  turbulence.  This  is  done  by  insuring  that  the 
modeled  evolution  equations  have  the  same  fixed  points  as  those  obtained  from  computer  and  labo¬ 
ratory  experiments  for  the  homogeneous  shear.  Results  of  computations  of  the  homogeneous  shear, 
with  and  without  rotation,  and  with  stabilizing  and  destabilizing  curvature,  are  shown.  Results 
are  consistently  better,  in  a  wide  class  of  flows  which  the  model  not  been  calibrated,  than  those 
obtained  with  other  nonlinear  models. 
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NA.S1-19480  while  the  author  was  in  residence  at  the  institute  for  Computer  Applications  in  Science  and  Engineering 
(1CASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681. 
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1.  Introduction 

Most  turbulence  models  are  devised  for  use  in  inertial  coordinate  systems.  Ad  hoc  changes  are 
then  made  to  reflect  the  unusual  effects  on  turbulence  of  swirl,  curved  streamlines  or  rotation  of 
the  coordinate  system:  one  makes  the  length  scale  depend  on  Richardson  number  or  adds  terms 
to  the  dissipation  equation  calibrated  to  observed  behavior.  This  approach  does  not  make  use  of 
the  mathematical  requirements  that  the  dependent  variables  and  their  evolution  equations  must 
satisfy  and  leads  to  models  that  perform  poorly  when  used  in  situations  substantially  different 
from  the  benchmark  flows  for  which  they  have  been  calibrated.  There  is  no  reason  a  second-order 
modeling  method  cannot  be  successfully  applied  to  a  high  Reynolds  number  rotating  turbulence: 
no  new  unknown  terms  appear  in  the  equations  and,  for  moderate  Rossby  number,  very  little  of 
the  phenomenology  on  which  the  technique  is  based  changes. 

The  effect  of  rotation  on  the  second-order  moments  is  felt  through  the  rapid-pressure-velocity  cor¬ 
relation  and  the  Coriolis  terms.  The  difficulty  with  the  equations,  as  modeled  presently,  can  be  seen 
when  the  equations  are  transformed  to  a  rotating  coordinate  system:  they  are  not  materially-frame- 
indifferent  in  the  two-dimensional  limit,  Speziale  (1981),  Hide  (1977).  This  phenomenon,  predicted 
by  the  Taylor- Proudman  theorem  requires  that  when  the  velocity  field  is  two-dimensionalized  by 
rapid  rotation,  without  components  along  the  axis  of  rotation,  the  equations  must  be  independent 
of  rotation.  The  problem  with  the  current  modeled  second-order  equations  results  from  the  inabil¬ 
ity  of  the  rapid-pressure  correlation  model  to  reflect  the  physics  embodied  in  the  Taylor- Proudman 
theorem.  This  is  another  form  of  realizability:  one  must  obtain,  for  a  bounded  flow,  in  the  limit  of 
rapid  rotation,  a  frame-indifferent  turbulence.  In  such  a  turbulence  the  rapid-pressure  correlation 
appearing  on  the  right  hand  side  of  the  second  order  moment  evolution  equations  equals  the  Cori¬ 
olis  terms  appearing  on  the  left  hand  side  leaving  the  equations  frame-indifferent.  A  rapid-pressure 
model  consistent  with  these  facts,  reflecting  more  of  the  information  contained  in  the  Navier-Stokes 
equations,  is  required. 

As  a  general  tool  is  being  developed  to  compute  a  wide  class  of  flows  for  which  there  may  not  be  any 
well  documented  benchmark  flows  with  which  to  calibrate  coefficients,  it  is  necessary  to  construct  a 
model  from  first  principles,  incorporating  more  of  the  physics.  In  the  present  new  representation  for 
the  rapid-pressure  strain  correlation  the  use  of  all  constants  that  do  not  come  from  first  principles 
have  been  minimized:  This  is  done  this  by  requiring  the  representation  to  have  the  proper  behavior 
in  five  different  limits:  1)  frame  invariant  with  respect  to  rotation  when  the  eigenvalue  of  the 
Reynolds  stress  tensor,  <  UiUj  >,  along  the  axis  of  rotation  vanishes,  2)  the  realizable  limit  in 
which  an  arbitrary  eigenvalue  of  the  Reynolds  stress  vanishes,  3)  the  joint-realizable  limit  in  which 
an  eigenvalue  of  the  tensor  <  99  ><  U{Uj  >  -  <  ><  On,  >  vanishes,  4)  the  isotropic  limit  in 

which  the  anisotropy  tensor,  6tJ  =  0,  and  5)  asymptotic  consistency  with  a  stationary  state  of  the 
turbulence  in  which  D/Dt  bij  =  0.  Of  these  five  principles  the  last  one  is  a  statement  based  on 
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experimental  observation  rather  than  mathematical  fact.  It  is  tacitly  asusmed  that  for  a  specific 
class  of  flows,  to  a  reasonably  suitable  approximation,  there  exists  an  equilibrium  state  to  which 
the  flow  relaxes  upon  the  removal  of  any  disturbing  forces. 

Assuring  the  proper  frame-invariance  of  the  modeled  second  moment  equations  is  done  by  requiring 
that  in  the  limit  of  a  turbulence  two-dimensional  in  planes  perpendicular  to  the  axes  of  rotation 
the  rapid-pressure  velocity  correlations  in  the  heat-flux  and  Reynolds  stress  equations  satisfy  the 
’’geostrophic”  constraints,  Ristorcelli  (1987)  and  Ristorcelli  and  Lumley  (1991b): 

tpqn^nXipqj  =  tqjn^n  ^  U{Uq  > 
tpqnQnX pqi  =  Vq&  ^  • 

The  term  geostrophic  is  borrowed  from  the  meteorological  literature  where  it  is  used  to  describe 
the  low  Rossby  number  balance  between  Coriolis  and  pressure  forces  in  the  evolution  equations  for 
large  scale  weather.  The  constraint,  however,  is  independent  of  how  the  flow  is  two-dimensionalized 
and  is,  therefore,  independent  of  Rossby  number.  The  requirements  of  the  Taylor- Proudman 
theorem  are  subsumed  by  the  requirements  of  the  principle  of  two-dimensional  frame  invariance 
(2DMFI)  first  put  forward,  in  the  context  of  turbulence  modeling,  by  Speziale  (1981).  A  flow  that 
is  two-dimensionalized  due  to  stable  stratification  or  magnetic  forces  is  also  frame-indifferent  in 
the  two-dimensional  two-component  limit  and  the  rapid- pressure  correlation  must  still  satisfy  the 
geostrophic  constraint  for  arbitrary  rotation. 

The  requirement  of  two-dimensional  material-frame-invariance  comes  from  first  principles:  it  is 
a  rigorous  limit  of  the  Navier-Stokes  equations  for  the  fluctuating  velocity,  Speziale  (1990).  The 
principle  of  2DMFI  seems  to  have  been  given  short  shrift.  Reynolds  (1989),  in  a  rapid  distortion 
theory  analysis  of  a  decaying  turbulence,  has  shown  that  rotation  has  a  distinct  effect  on  the  sec¬ 
ond  invariant,  II,  of  the  anisotropy  tensor,  bij,  and  goes  on  further  to  show  that  the  equations  for 
the  invariants  using  the  current  rapid-pressure  models  are  independent  of  rotation.  He  concludes 
that  such  models,  because  they  don’t  show  an  explicit  dependence  on  the  rotation  in  the  evolu¬ 
tion  equations  for  the  invariants,  are  inadequate  for  rotating  flows.  A  reference  to  the  evolution 
equations  for  the  invariants,  fc;  example  Speziale  et  al.  (1991),  shows  that  both  the  production 
and  dissipation,  which  are  functions  of  the  rotation  through  their  dependence  on  the  Reynolds 
stresses,  appear  in  the  evolution  equations  for  the  invariants.  The  equations,  except  for  the  case  of 
a  decaying  turbulence  in  which  there  are  no  production  terms,  are  not  independent  of  rotation. 

Shih  and  Shabbir  (1990),  following  Reynolds  argument  decide  that  since  the  rotation  terms  do  not 
explicitly  appear  in  the  evolution  equation  for  the  invariants  (which  is  true  for  any  rapid-pressure 
model  that  is  a  function  of  symmetric  tensors)  that  the  use  of  a  properly  frame-indifferent  model  will 
not  produce  the  proper  turbulence  behavior.  This  argument  is  specious:  the  same  argument  applied 
to  the  trace  of  the  Reynolds  stress,  the  kinetic  energy  of  the  turbulence,  q2  =<  uvup  >,  one  would 
abandon  refinement  of  rapid-pressure  models  because  they  do  not  affect  the  energy  of  the  turbulence 
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-  the  evolution  equation  for  q 2  being  independent  of  the  rapid- pressure  modeling.  However,  different 
rapid-pressure  models  produce  different  Reynolds  stresses  and  different  turbulence  energies  in  the 
same  flows.  Even  if  the  models  fails  to  account  for  some  predictions  of  RDT  for  a  restricted  class 
of  flows  neglecting  results  that  are  rigorous  consequences  of  the  Navier-Stokes  equations,  Speziale 
(1990),  are  not  justified.  (It  should  be  kept  in  mind  that  the  RDT  limit  violates  various  equilibrium 
assumptions  made  in  the  mathematical  development  of  the  second-order  modeling  method).  The 
2DMFI  constraint  is  relevant  to  flows  other  than  rotating  flows:  any  flows  in  which  the  largest 
scales  of  the  motion  are  two-dimensionalized  by  strong  body  forces  will  come  arbitrarily  close  to 
being  rotationally  frame  invariant.  A  dramatic  example  of  this  might  be  the  collapse  of  the  wake 
in  a  stably  stratified  environment.  Shih  and  Shabbir  (1990)  also  cite  Reynolds  (1987)  as  having 
shown  that  the  2DMFI  constraint  does  not  provide  any  new  information.  Here  it  is  found  that 
the  geostrophic  constraint  produces  six  additional  linearly  independent  constraint  equations  for  the 
rapid-pressure  model. 

The  present  article  is  a  description  of  the  derivation  and  validation  of  a  rapid-pressure  model 
that  satisfies  the  principle  of  material-frame-indifference  in  the  two-dimensional  limit.  The  present 
derivation  of  the  2DMFI  model  has  been  given  earlier,  in  departmental  reports  in  Ristorcelli  (1987, 
1991)  or  Ristorcelli  and  Lumley  (1991b).  In  those  previous  developments  of  the  2DMFI  model 
several  free  parameters  were,  for  simplicity,  set  to  zero.  In  this  paper  the  stationary  points  of  the 
homogeneous  shear,  as  reflected  in  the  structural  equilibrium  assumption,  D/Dt  bl}  =  0  where 
b{j  =<  U{Uj  >  /  <  upup  >  -  1/3 6{j,  are  used  to  set  the  unknown  free  parameters  occurring  in  the 
model.  The  algebraic  equations  describing  the  stationary  states  of  the  homogeneous  shear  along 
with  the  experimentally  determined  values  of  the  stationary  state  are  used  as  constraints  to  set 
the  free  parameters  in  the  model.  This  insures  that  the  fixed  points  of  the  modeled  differential 
equations  are  the  same  as  the  fixed  points  of  a  particular  physical  flow.  This  calibration  is  done 
without  sacrificing  any  of  the  mathematical  principles  built  into  the  model. 

The  next  section  of  this  article  defines  the  problem.  A  third  section  presents  a  derivation  of  the 
tensor  polynomial  model  for  the  rapid-pressure  correlation.  A  discussion  of  the  constraints  as  well 
as  a  derivation  of  the  new  geostrophic  constraint  is  given.  A  new  ansatz  for  the  rapid-pressure 
correlation  produces  a  model  valid  away  from  both  the  geostrophic  (two-dimensional  in  planes 
perpendicular  to  the  axis  of  rotation)  and  realizability  limits  (two-dimensional  in  that  one  of  the 
eigenvalues  of  the  Reynolds  stress  tensor  dissappears).  The  issues  of  realizability  in  the  light  of 
some  recent  work,  Speziale  (1993),  are  discussed.  A  subsequent  section  uses  results  from  matrix 
algebra  to  collapse  the  model  to  a  structure  similar  to  other  rapid-pressure  models  with  some 
interesting  differences.  The  model  is  seen  to  have  the  same  tensor  bases  as  the  FLT  model  and  in 
the  limit  of  a  planar  mean  flow  it  has  the  same  tensor  bases  as  the  quasi-linear  SSG  model.  The 
two  subsequent  sections  show  computations  done  with  the  model  and  compare  the  results  to  other 
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models.  Following  this  a  general  discussion  of  modeling  issues,  in  the  context  of  realizability  and 
the  use  of  stationary  states  as  the  only  self-consistent  asymptotic  limit  for  calibrating  turbulence 
models  is  given. 

2.  The  rapid-pressure  correlation 

In  an  incompressible  turbulence  in  a  rotating  coordinate  system,  with  buoyancy  effects  in  the 
Boussinesq  approximation,  the  Reynolds  stress  equations  have  the  following  form: 

DfDt  <  tijUj  >  +2 (tifcp  <  upUj  >  +(jkP  <  UpUi  >)f l*Ro-1  =  +  <  Oui  >  /3j+  <  Otij  >  /3, 

U  J lip  f/, ,p  *f  ^  U, Up  ^  U^Ujtlp 

Pij  u^  ^  "f*  ^  p,i  Hj  Re  ^  UjUj  2 Re  ^  ujyp  ^ 

D/Dt  <  0u,  >  +2<pifcfljt  <  0up  >  Ho-1  =  -[<  0Uj  >  f/tiJ  +  <  UiUj  >  T,j]+  <  00  >  ft 

-  <  OuiUj  >j  -  <  p  ,0  > 

+  Re_1(l  +  Pr_1)(<  0U{  >jj  -2  <  9,jUij  >  . 

The  velocity  has  been  normalized  by  a  characteristic  velocity  uc  and  the  Rossby  number  is  Ro  = 
uJSlRc  where  Rc  is  a  length  scale  and  Q.  the  rotation  rate  of  the  frame  of  reference.  The  gravity 
and  rotation  vectors  are  aligned  with  the  3  axis.  Our  concern  is  with  the  pressure- velocity  and 
-temperature  correlations,  <  p,j  u,  >  and  <  p„  0  >.  An  equation  for  the  pressure  fluctuations 
comes  from  the  divergence  of  the  Navier-Stokes  equations  for  the  fluctuating  velocity 

uiit  —  <  U{Uj  >,j  +2 C,|tpfl)cUpfto  1  =  ~P>i  4"^ft  +  Re  1  U i,jj 

which  produces  a  Poisson  equation  for  fluctuating  pressure.  The  standard  linear  decomposition 
recognizes  three  terms 

~Pr,ii=  2[Ui,p+tpiktokRo~1]up,t 
~P  »i  =  Ui  ij  uj  ii  “  ^  ui  ij  uj  >«  > 

P,ii  =  PiO,j 

where  pT,p3,pb  are  respectively  the  rapid-pressure,  the  slow  or  return  to  isotropy  pressure,  and  the 
buoyancy-pressure.  The  effects  of  rotation  are  felt  through  the  rapid-pressure,  pr .  Solution  of  the 
Poisson  equation  for  the  rapid-pressure  is  by  application  of  Green’s  theorem 

4>{x)  =  -  (47T)-1  /</>(x'),jjdx7(x  -  x'). 

It  is  the  moments  of  the  solution  that  are  required  to  close  the  second-order  equations.  For  a 
homogeneous  mean  field,  more  than  an  integral  scale  away  from  any  surfaces,  a  staightforward 
interchange  of  the  order  integration  and  averaging  produces: 

K  P  ,j  Ui  >  +  <  p  ,t  Uj  >=  —  2[f/,j,p  -\-(pqk^ik  Ro  ^  ]  [Xipqj  +  Xjpqj] 

<  Pr,i9  >=  -2[Uq,p  +(pqktokRo~l)Xpiq 

where 

Xpiq  =  (47r)-1  /  <  0(x)Up(x')  >„y  dx’/(x  -  x1) 

Xipqj  =  (4*)-1  /  <  Uj(x)Up(x')  > ,j'q>  dx 7(x  -  x'). 
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The  construction  of  a  tensor  polynomial  model  for  the  volume  integrals  of  the  two- point  correlations, 
Xijki  and  Xijk,  is  the  subject  of  this  article.  The  rapid-pressure  covariance  integral  appearing  in 
the  heat  flux  equations,  Xijk ,  is  treated  as  the  Reynolds  stresses  and  the  heat  fluxes  are  linked 
through  Cauchy-Schwarz  inequalities,  <  BB  ><  UiUj  >  -  <  Bu,  ><  Buj  >  >  0  and  their  modeling 
cannot  be  done  independently.  The  rapid- pressure  term  appearing  in  the  Reynolds  stress  equations 
necessary  to  treat  the  mechanical  turbulence  problem  is  the  focus  of  this  article. 

3.  Obtaining  a  representation  for  the  rapid- pressure  correlation 

A  constitutive  relation  in  which  the  integral  of  the  two-point  covariance  is  parameterized  by  a  local 
function  of  the  anisotropy  tensor  and  heat-flux  vector  is  proposed.  The  most  general  forms  of  such 
relationships  are  the  following  tensor  polynomials 

Xijki/  <  upup  >=  A.\ SijSki  4  A2(6ikSji  4  SuSjk) 

4  AsSijbki  4  A+bijSki  4  A^(bikS3i  4  buSjk  4  Sikbji  4-  Subjk) 

+  A6St]blt  4  Ajb^Ski  4  Ag(b?k6ji  4  b^Sjk  4  6, *6^  4  f>ubjk) 

4  Agbjjbki  4  Aio(bikbji  4  bubjk) 

4  Anbijbh  4  6*/  4  A\3{b2kbji  4  blffck  4  4  bub2jk) 

4  Aub’jbl  4  AntfS  4  b2b%) 


Xpkj  =  Di  <  0up  >  Skj  4  D2(<  6uk  >  Spj 4  <  Bu:  >  6pk) 

4  D3  <  Bup  >  bkj  4  Di(<  Buk  >  bPJ ;4  <  Buj  >  bpk) 

4  Ds  <  Bup  >  b2kj  4  D6(<  Buk  >  bj}+  <  Buj  >  b^) 

4  [D7bqp6kj  4  Dg(bqk6pj  4  <  ®uq  > 

4  [D9bqpbkj  4  Dio(bqkbpj  4  bqjbpk))  <  Buq  > 

+  [Du bqpb\j  4  Dl2(bqkblj  4  bqjblk)}  <  Buq  > 

4  [Di3b2bkj  4  D\4(b2k6pj  4  b2j6pk)]  <  Buq  > 

4  [Dnbqpbkj  4  D\g{b,bpj  4  b  jbpk)]  <  Buq  > 

4  [Dnb2pbkj  4  Dxg{b)kb2V]  4  %%)]  <  6uq  > 

where  btJ  is  the  anisotropy  tensor  6^  =<  «,Uj  >  /  <  uquq  >  -1/3 6ij  and  <  Bui  >  is  the  turbulent 
heat-flux.  Following  Pope’s  linearity  principle  only  terms  linear  in  the  heat-flux  are  kept.  The  A, 
and  Di  are  functions  of  the  invariants  of  6,j  and  <  Bui  >. 

Parameterizing  the  integral  of  a  two-point  correlation  in  terms  of  the  local  anisotropy  tensor  is  a 
substantial  simplification  requiring  consideration.  For  a  turbulence  with  short  term  memory  and 
limited  awareness  Lumley  (1970)  has  discussed  the  conditions  under  which  such  a  constitutive 
relation  is  tenable  and  carried  out  a  similar  expansion  procedure,  Lumley  (1967),  indicating  how 
the  truncation  errors  scale.  ^From  one  point  of  view,  the  constitutive  relation  proposed  can  be  seen 
as  the  first  term  in  a  functional  Taylor  series  expansion  for  the  rapid-pressure  correlation.  As  the 
correlation  decays  with  distance  the  primary  contribution  to  the  integral  will  come  from  regions 
within  an  integral  length  scale  of  the  local  position  -  thus,  in  a  homogeneous  turbulence,  the  first 
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term  will  constitute  a  good  approximation.  Retaining  higher  order  terms  of  the  functional  Taylor 
series  expansion  which  involves  spatial  and  temporal  derivatives  of  bt}  substantially  complicates 
the  problem.  It  is  expected  that  the  retention  of  only  the  first  order  term  captures  enough  of  the 
physics  to  allow  prediction  suitable  for  engineering  purposes. 

In  applying  the  above  constitutive  relation  in  non-inertial  systems  there  are  several  phenomena, 
peculiar  to  rotating  flows,  that  may  be  in  conflict  with  the  assumptions  of  the  turbulence  closure. 
In  particular  the  inertial  wave  field  associated  with  the  rotation  may  interfere  with:  1)  the  energy 
cascade  from  the  large  to  the  small  scales  of  the  flow,  2)  the  universal  equilibrium  assumed  for 
the  small  scales  of  the  flow,  and  3)  the  assumed  steadiness  of  the  mean  flow.  These  issues  require 
consideration. 

In  a  mean  flow  rotating  with  speed  ft  there  will  be  inertial  oscillations  with  frequency  less  than  or 
equal  to  2ft.  For  a  quasi-steady  assumption  to  be  valid,  changes  in  the  mean  flow  must  be  slow  with 
respect  to  the  turbulence’s  ability  to  adjust  to  the  imposed  changes.  This  requires  that  2ft  <  e/k 
and  becomes  a  lower  bound  for  a  Rossby  number,  Ro  =  e/2  ilk  >  1,  for  which  an  equilibrium  theory 
is  appropriate. 

Similar  to  a  stably  stratified  density  field  in  which  stratification  inhibits  particle  motions  in  the 
vertical  direction,  rotation  inhibits  transverse  displacements  of  fluid  particles.  Similar  to  the  radius 
of  gyravion  of  a  charged  particle  in  a  magnetic  field,  the  lateral  displacement  of  a  turbulent  fluid 
element  in  a  rotating  frame  can  be  characterized  by  alength  scale  (g2/3)1/,2/2ft  where  q2  =<  UjU,  >. 
This  length  scale  must  be  larger  than  the  turbulence  length  scale  l  ~  (q2/ 3)3/,2/e  to  insure  that  the 
transverse  confinement  by  the  Coriolis  forces  does  not  affect  the  fluctuating  field.  This  produces  a 
similar  bound  on  the  turbulent  Rossby  number  Ro  —  e/2 ftfc  >  3/2. 

The  phase  coherence  necessary  for  the  cascade  of  energy  to  the  smaller  scales  of  the  motion  will 
be  interfered  with  if  the  production  scales  of  the  motion  k£  ~  1,  where  n  is  the  wavenumber,  are 
subject  to  an  inertial  wave  field.  A  different  Rossby  number  can  be  defined  as  a  ratio  of  the  vorticity 
of  the  production  scales  of  the  motion  to  the  background  vorticity  Rot  =  (g2/3)J/2/2ft/  =  3e/2 Slq2 
using  e  =  (q2/3)1/,2/f.  A  spectral  Rossby  number  can  also  be  defined  as  Ro(k )  =  u(k)/2SU(k)  = 
(kE(k))1/2/2SI(2w /k)  which  using  the  inertial  range  scaling  E(k )  =  ae2/3*'5/3  and  e  =  (q2/3)1|/2/f 
becomes  Ro(k)  ~  0.2 (Ki)2^3Rot.  The  effects  of  rotation  decrease  as  the  wave  number  increases.  For 
the  inertial  oscillations  associated  with  the  rotation  not  to  interfere  with  the  cascade  mechanism 
Ro(k)  >  1  for  ni  ~  1  is  required.  Thus  for  Rot  >  5  the  usual  parameterization  of  the  spectral 
cascade  rate,  e,  in  terms  of  the  energy  containing  scales  of  the  motion  is  appropriate.  For  Rot  <  5  the 
current  dissipation  equation  requires  modification.  How  the  dissipation  equation  is  to  be  changed 
to  account  for  the  effects  of  rotation  on  the  cascade  rate  is  an  unresolved  issue  that  is  the  topic 
of  current  research.  It  is,  however,  clear  that  the  assumption  of  the  small  scale  equilibrium  with 
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the  large  scales  of  the  motion  is  valid  in  most  high  Reynolds  number  rotating  flows  of  interest: 
the  Rossby  number  of  the  dissipation  scales  of  the  motion  is  Roe  =  Re1/*  Rot  and  therefore  the 
dissipation  scales  are  insensitive  to  the  effects  of  rotation. 

Having  dealt  with  some  of  the  phenomenological  issues,  the  mathematical  constraints  required  of 
the  pressure-velocity  correlation  are  now  investigated. 

There  are  several  physical  inequalities  and  mathematical  identities  the  X  must  satisfy.  These  prin¬ 
ciples  are  used  to  obtain  a  set  of  algebraic  constraint  equations  for  the  Ai  and  Z),.  For  an  arbitrary 
three-dimensional  turbulence  the  tensor  polynomials  must  satisfy  the  symmetry  constraints, 

Xijkl  —  Xijlk 

Xijkl  —  Xjikl 

Xijk  =  Xikj- 

These  symmetry  constraints  are  built  in  to  the  assumed  form  of  the  tensor  polynomials.  For  an 
arbitrary  three-dimensional  turbulence  the  tensor  polynomials  must  also  satisfy  the  constraints  of 
normalization,  continuity: 

Xtjkk  =  <  U««j  >  Xikk  =  <  ««0  > 

Xijjk  =  0  Xiji  =  0. 

Note  that  a  contraction  of  the  integral  of  a  two-point  statistic  is  a  local  one-point  statistic. 

The  tensors  <  u,Uj  >  and  <  60  ><  U{Uj  >  -  <  0Ui  ><  0uj  >  are  positive  semi-definite.  This 
reflects  the  fact  that  the  energy  of  the  turbulence  is  always  positive  and  that  the  magnitude  of  the 
correlation  coefficients  between  the  various  components  of  the  tensors  be  bounded  by  one.  These 
facts  lead  to  the  ’’realizability”  and  "joint-realizability”  constraints  which  specify  the  behavior  of 
the  correlations  when  specific  limit  states  are  approached.  The  relevant  portion  of  the  Reynolds 
stress  transport  equations,  in  principal  axes,  requires  that 

D/Dt  <  uaua  >  ~  [Up,i  +(ipk^kRo~1]Xiapa  -*  0  as  <  uaua  >-►  0 

in  order  to  satisfy  realizability.  The  rate  of  change,  due  to  the  rapid-pressure  correlation,  of  the 
eigenvalue  <  uaua  >  is  required  to  vanish  as  the  limit  state  is  approached.  This  insures  that  the 
rapid-pressure  correlation  model  does  not  cause  the  solution  to  go  into  the  unrealizable  region  in 
which  <  uaua  >  is  negative.  This  realizability  limit  is  rephrased  in  terms  of  the  determinant  of 
the  Reynolds  stress:  F  =  (R?--  -  3RjjRF  +  2R^)/6  where  R,j  =<  >  /  <  uPup  >  which 

can  be  written  in  terms  of  the  invariants  of  the  anisotropy  tensor  as  F  =  1  +  91 1  +  271 1 1  where 
II  =  —\/2bijbij  =  -1/2  <  b 2  >,  III  =  l/3bipbpjbji  =  1/3  <  63  >.  The  determinant  F  varies 
between  zero  and  one;  F  =  1  corresponds  to  an  isotropic  turbulence  and  F  =  0  corresponds  to  the 
realizable  limit. 
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Similar  reasoning  applied  to  the  mixed  tensor  involving  the  Reynolds  stress,  the  heat  flux  and  the 
variance  of  the  temperature  fluctuations,  produces  the  "joint- realizability”  constraint 

D / Dt  Daa  ~  [^p,i  -bCipfcHfc Ro  '][<  99  >  Xjapa~  <  9ua  >  Xiap]  ♦  0  as  Daa  •  0 

which  couples  the  rapid-pressure  correlations  appearing  in  the  heat  flux  and  the  Reynolds  stress 
equations.  A  similar  determinant  function  Fj  is  defined  with  the  normalized  D{},  for  which  0  < 
Ffi  <  \.  Joint- realizability  reflects  the  requirement  that  the  magnitude  of  the  correlation  coefficients 
be  bounded  by  one:  the  Reynolds  stress  and  the  heat-flux  take  on  values  "jointly”  such  that  the 
time  rate  of  change  of  Daa  goes  to  zero  as  Daa  goes  to  zero.  Additional  detail  regarding  the 
application  of  realizability  constraints  can  be  found  in  Shih  and  Lumley  (1985). 

Note  that  no  assumptions  regarding  the  higher  order  derivatives  of  the  eigenvalue  have  been  made. 
The  strong  form  of  realizability,  in  which  D2/Dt2  <  uaua  >  >  0  is  required  at  the  realizability 
limit  in  order  to  allow  the  turbulence  to  leave  the  realizable  state  is  not  invoked.  Such  an  agency  is 
already  present  in  the  slow  terms  and  it  is  not  necessary  to  force  the  rapid  terms  to  be  responsible 
for  such  behavior  which,  as  will  be  seen,  is  inconsiitent  with  the  a  small  parameter  expansion 
of  the  rapid-pressure  representation  around  the  realizability  limit.  Instead  a  weak  realizability 
constraint,  as  specified  by  Speciale  et  al.  (1993),  Pope  (1983),  which  does  not  allow  the  solution 
to  attain  the  realizable  limit  in  finite  time  is  invoked.  This  is  done  by  requiring  that  the  rapid- 
pressure  correlation  vanish  more  rapidly  than  the  slow-pressure  correlation  model.  This  avoids 
any  assumptions  regarding  the  behavior  of  the  second  derivative  which  are  required  for  the  flow  to 
leave  the  realizable  state  which  is  accessible  in  finite  time  in  models  using  the  strong  form  of  the 
realizability  constraint.  Moreover  recent  work  by  Speziale  et  al.  (1993)  indicates  that  the  present 
hierarchy  of  second-order  models  is  inconsistent  with  the  strong  form  of  the  realizability  constraint. 
The  rate  of  rotation  of  the  eigenvalues  coming  from  the  second  derivative  is  a  sink  term  of  a  form 
that  that  cannot  be  balanced  by  the  present  models.  In  setting  the  portion  of  D/ Dt  <  uaua  > 
due  to  the  rapid-pressure  correlation  to  zero  while  choosing  a  return  term  model  that  precludes 
accessibility  of  the  realizable  limit  state  the  weak  form  of  realizability  is  satisfied  by  the  sum  of  the 
modeled  terms  on  t  he  right  hand  side  of  the  transport  equations.  The  issue  raised  by  Speziale  et 
al.  (1993)  does  not  impact  on  the  present  model,  as  all  the  higher  order  derivatives  vanish,  and  the 
realizable  limit  is  not  attainable  in  finite  time  because  the  relative  rate  of  dissappearance  of 
slow-pressure  terms  with  respect  to  the  rapid-pressure  terms.  It  may,  however,  require  a  rethinking 
of  the  modeling  principles  so  that  the  intrinsic  negativity  of  the  second  derivative  can  be  properly 
balanced  if  the  realizable  limit  is  to  be  considered  accessible. 

The  last  of  the  constraints  to  be  applied  to  specify  the  constitutive  relation,  the  geostrophic  con¬ 
straint,  is  now  discussed.  The  derivation  of  the  constraint,  as  it  has  not  been  given  elsewhere  in 
the  published  literature,  Ristorcelli  (1991),  is  given.  For  a  turbulence  two-dimensional  in  planes 
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perpendicular  to  the  ax’s  of  rotation  Hide(1977)  and  Speziale  (1985,  1990)  have  shown  that  the 
equations  are  materially- frame- indifferent.  This  is  a  direct  consequence  of  the  Taylor- Proudman 
theorem.  For  a  modeled  set  of  equations  to  be  consistent  with  the  2DMFI  principle  the  rapid- 
pressure  correlations  must  satisfy  the  "geostrophic”  constraints, 

tpqnQn%ipqj  =  tqjn^n  <  UjUq  > 
tpqn^nXpqi  =  (, jin  fin  ^  Uq6  >  . 

Consider  the  portion  of  the  rapid-pressure  correlation  associated  with  the  rotation 

tpqkftkRo~l  Xipqj  =  (pqknk(4*Ro)~l  /  <  U;(x)tip(x')  <fx'/(x  -  x'). 

A  velocity  field,  two-dimensional  in  planes  perpendicular  to  the  axis  of  rotation,  has  the  representa¬ 
tion  Up  =  tpqkSlk4>,q-  Inserting  the  expression  for  the  velocity  field  into  the  integral  and  contracting 
produces,  in  the  integrand,  the  Laplacian  of  the  streamfunction,  i>,p=  tqpkSlkuq  which  reduces  the 
volume  integral  of  a  two-point  statistic  to  a  local  one-point  statistic, 

tpqk^k  Ro  X ipqj  —  Ro  C  Uj0,j  ^  —  Ro  ^pjkQk  ^  U,Up  ^ 

upon  application  of  Green’s  theorem.  The  geostrophic  constraint  €pgnfl„.Xjpgj  =  €gj„fl„  <  utuq  > 
then  follows.  The  satisfaction  of  the  geostrophic  constraint  means  that,  when  the  turbulence  is 
two-dimensional  in  planes  perpendicular  to  the  axis  of  rotation,  the  portion  of  the  rapid-pressure 
correlation  associated  with  the  rotation  is  equal  to  the  Coriolis  terms  effectively  removing  any 
dependence  on  the  rotation  rate  from  the  equations.  Such  a  Coriolis-pressure  force  balance  is 
found  in  the  large  scale  atmosphere  were  it  describes  the  well-known  geostrophic  wind.  A  point  on 
nomenclature:  a  geostrophic  turbulence  is  one  that  is  two-dimensional  because  its  Rossby  number 
is  small.  This  latter  qualification  distinguishes  it  from  a  two-dimensional  turbulence  (that  is  made 
so  by  some  other  means)  for  arbitrary  Rossby  number.  Both,  however,  are  frame-indifferent. 

The  application  of  the  five  sets  of  constraints  -  normalization,  continuity,  realizability,  joint- 
realizability  and  geostrophy  -  produce  thirty-six  linear  algebraic  equations  (several  of  which  are 
redundant)  for  the  thirty-three  unknown  coefficients  Ai,(i  =  1,15)  and  D„ (t  =  1,18)  appear¬ 
ing  in  the  tensor  polynomials.  The  equations  are  of  the  general  form  A,y(//,  III)xj  =  bi  where 
xi  —  [Ai,...,  A15,  Djg]  and  where  II  and  III  are  the  invariants  of  the  anisotropy  tensor, 

//  =  -l/26fj6,j  =  —1/2  <  b 2  >,  III  —  l/3bivbvjbji  =  1/3  <  b3  >.  Using  the  definition 
F  =  1  +  9//  +  27/ 1 1  the  general  form  of  the  equations  can  be  rewritten  as  A,  _,(//,  F)xj  =  6.. 
Note  that  F  and  II  appear  linearly  in  the  constraint  equations. 

The  ansatz 

Xi,u=  Ku  +  FX? 

Xiik=  xflt  +  Fxfjk 

is  used  to  extract  more  information  from  the  constraint  equations.  Here  X°  satisfies  the  set  of 
constraint  equations  A,y(//,0)x®  =  6,  obtained  by  application  of  all  five  sets  of  constraints  which 
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are  simultaneously  valid  only  when  F  =  0  (or  equivalently  when  III  =  -(//  +  l/9)/3).  While 
XF  is  obtained  from  the  reduced  set  of  equations  A,y(//,  F) xf  =  6®  which  satisfy  the  three  sets  of 
constraints  of  normalization,  continuity  and  joint-realizability  and  where  6®  is  a  known  function  of 
the  X®  solution.  As  there  are  more  unknowns  than  equations  there  are  free  parameters  which  will 
be  used  shortly  to  insure  that  the  model  is  asymptotically  consistent  with  an  equilibrium  state. 
The  solution,  in  which  all  free  parameters  are  set  to  zero,  and  which  satisfy  all  the  mathematical 
constraints,  will  be  called  the  basic  model  and  has  been  given  in  Ristorcelli  (1987,1991)  and  Ris- 
torcelli  and  Luniley  (1991b).  These  coefficients  are  given  in  Appendix  3.  It  is  to  this  basic  model 
that  additional  terms  are  added  in  order  to  insure  that  the  model  is  asymptotically  consistent  with 
a  known  equilibrium  state  of  a  particular  turbulent  flow  field. 

In  reference  to  the  coefficients,  A,,  of  the  basic  model  and  defined  in  Appendix  3,  several  comments 
are  appropriate.  The  A,  are  nonlinear  functions  of  the  invariants  of  btJ:  they  are  ratios  of  polyno¬ 
mials  of  the  invariants.  The  realizability  limit  is  attained  along  the  line  III  =  -(II  +  1  /9)/3  on 
which  F  —  0.  At  the  isotropic  limit  II  =  III  =  0  and  F  =  1,  and  thus  A\  -f  A2  =  1/10  and  the  well 
known  exact  result  for  isotropic  turbulence  is  obtained.  Note  that  the  ”  off-realizability”  correction 
FXf  is  necessary  to  obtain  this  limit.  The  satisfaction  of  the  isotropic  limit  is  discovered  to  be 
a  consequence  of  the  constraints  used  to  create  the  2DMFI  model  -  it  is  not  a  constraint  that  has 
been  enforced  to  obtain  the  model  but  is  satisfied  naturally  by  the  ansatz.  This  seems  to  vindicate 
the  present  procedure. 

Note  that  the  coefficients  appear  to  be  singular  at  the  one-dimensional  limit  when  1  +  3//  =  0. 
This  is  actually  not  the  case  as  the  singularities  arising  from  the  individual  terms  in  the  rapid 
pressure  representation  annihilate  each  other  when  summed.  This  has  beeen  verified  analytically 
near  the  one  dimensional  limit,  as  II  — >  -1/3,  in  principle  axes  of  the  Reynolds  stress  tensor  using 
a  perturbation  expansion. 

It  should  be  emphasized  that  the  form  of  the  rapid-pressure  correlations  was  chosen  so  that  the 
coefficients  are  valid  for  all  states  of  the  turbulence  -  they  are  not  fixed  to  their  values  at  the  real¬ 
izable  limit.  The  resulting  model  therefore  satisfies  all  mathematical  constraints  for  any  arbitrary 
Reynolds  stress  -  not  just  at  the  realizability  or  geostrophic  limit  states.  The  second  term,  FXf, 
in  the  expression  for  X,  represents  the  ”off-realizability  correction”,  thus  the  coefficients,  A,  and 
Dt,  are  functions  dependent  on  the  state  of  the  turbulence  as  parameterized  by  the  invariants  and 
only  attain  their  realizable  limit  values  at  F  =  0. 

4.  A  compact  representation  of  the  2DMFI  rapid-pressure  model 

The  rapid-pressure  correlation  models  are  usually  written  in  the  form  they  appear  in  the  Reynolds 
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stress  equations  with: 


IE  =  2[lf„p+f  pqkftkRo^jXpi,, 

n«j  =  2 [Uq,p +tpqk(lkRo  +  Xjpqi}- 

which  can  be  rewritten  in  terms  of  the  strain  and  rotation  tensors  as 


Uij  —  2[5?p  +  Wqp)[Xipq  j  +  Xjpqi] 

where  Sqp  —  1/2 {Uq,p  +  Up,q)  and  —  [l/2(f/?iP  —  Up<q)  +  (pqif£lkRo~l]  are  the  usual  mean  strain 
rate  and  the  total  or  intrinsic  rotation  rate.  The  second-order  equations  are  then  rewritten  as 

D/Dt  <  U{Uj  >  +2 (ikp  <  Uptij  >  SlRo~x  +  2fjfcp  <  uptfj  >  HkRo~ 1  =  -I-  +••• 

d  /  t)t  <  Out  >  +2(pikSik  <  >  Ro~ 1  = + nr + nr  +— 

where  the  terms  omitted  have  already  been  given.  Taking  the  contraction  of  the  fourth-order  tensor 
on  the  mean  velocity  gradients  to  obtain  the  form  used  in  the  Reynolds  stress  equations  produces, 
with  q2  =<  Up  Up  >, 

m,  /2<?2  =  2[Ai  +  A2}Sii 

+  [(A3  +  A\  +  2As)(bipSpj  +  bjpSpi)  +  2AsS{j  <  bS  >] 

+  [(As  -  A4)(bipWpj  +  bjpWpi)) 

+  {(A*  +  A7  +  2A8)(bfpSpj  +  b%Spi )  +  2 AtSij  <  b2S  >] 

+  [(A,  -  A,)(b>pWpi  +  b)twt,)\ 

+  2[(A9  +  Aio)6ip5M6,j  +  Aiobjj  <  bS  >] 

+  [(An  +  A, 2  +  2Aa)(birS„b2,j  +  bipS„b2,)  +  2 A,j(i?.  <  bS  >  +b„  <  b‘S  »] 

+  [(An  -  A,2)(6jp W„bl  +  b„W^)] 

Here  the  angle  brackets  represent  the  trace  of  the  indicated  quantity:  eg.  <  bS  >=  bX]  StJ  and 
that  II  =  -1/2  <  b 2  >,  III  =  1/3  <  b3  >.  Note  that  the  n^  has  zero  trace  because  of 
the  continuity  constraint,  Xijjk  =  0,  requires  A3  +  A4  +  5A5  -  //(An  +  A\2  +  4A13)  =  0  and 
A6  +  A7  +  5Ag  +  A9  +  A10  =  0  It  is  possible  to  rewrite  the  higher  order  tensor  bases  in  terms  of 
the  lower  order  terms  substantially  simplifying  the  form  of  the  model.  The  generalized  Cayley- 
Hamilton  theorem  is  used  to  rewrite  the  expression  in  an  irreducible  tensor  basis.  Using  the  matrix 
notation 

bSb  =  -[b2S  +  Sb2]+  <  bS  >  b  +  1/2  <  b2  >  S+  <  b2S  >  1 
bSb2  +  b2Sb  =  -1/3  <  b3  >  S+  <  b2S  >  b+  <  bS  >  b2 

the  2DMFI  rapid-pressure  correlation  can  be  written  more  compactly  as 

=  [B3+  <b2>  BS+  <b3>  BflSij 

+  B4{bipSpj  +  bjpSpi  -  2/3  <bS>  Si,}  +  B'l'  <bS>  [6?  +  2///3dti] 

+  Bs[bipWpj  +  b^Wpi]  +  [fl6  <  bS  >  +B"  <  b2S  >]btJ 
+  B7[b2pSPJ  +  b%Spt  -  2/3  <  b2S  >  6tJ] 

+  B8[b2ipWpJ  +  b)pWpi]  +  B9[bipWqpb2i  +  b]pWqpb2x] 
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where  the  values  of  the  B,  in  terms  of  the  A{  are  given  in  Appendix  1 .  The  above  representation 
of  the  rapid-pressure  correlation  will  be  called  the  basic  or  uncalibrated  model.  The  coefficients, 
B,,  appearing  in  the  base  2DMFI  model  above  come  from  first  principles:  they  do  not  result  from 
any  numerical  optimization  with  experimental  or  numerical  data.  Comparisons  to  the  FLT  model, 
Fu  et  al.  (1987),  shows  that  the  two  models  have  the  same  tensor  structure  containing  the  same 
generators.  The  only  difference  is  in  the  terms  proportional  to  [6^]  and  [6?  4-  211  /36,,\,  usually 
identified  with  the  slow-pressure’s  contribution  to  the  pressure-strain  correlation  and  calibrated 
accordingly,  now  reflect  a  contribution  from  the  rapid-pressure.  Contraction  of  the  irreducible  form 
of  X, jki  on  the  mean  velocity  gradients  has  produced  the  bases  [6^]  and  [6^  +  211  /36t,}  in  which 
the  coefficients  are  functions  of  the  invariants  <  65  >  and  <  b2S  >  in  addition  to  the  dependence 
on  <  6 2  >  and  <  63  >  appearing  in  the  A{.  The  two  tensor  bases,  usually  associated  with  the 
slow-pressure  correlation,  arise  as  a  consequence  of  starting  with  the  two-point  volume  integral 
and  represent  a  contribution  to  the  pressure-strain  correlation  whose  structure  is  identical  to  the 
slow-pressure  models  but  whose  genesis  is  in  the  rapid-pressure. 

Speziale  et  al.  (1991)  have  also  written  a  general  form  for  the  pressure-strain  covariance.  It  is  linear 
in  the  mean  velocity  gradients  and  nonlinear  in  the  anisotropy  tensor.  Their  expression  contains 
the  same  generators  as  the  present  model  except  for  the  cubic  term  [bWb2  -  b2Wb].  Speziale 
et  al.  (1991)  use  results  from  rational  mechanics  (cf.  Smith  (1971))  to  expand  in  a  functional 
basis.  The  present  strategy  uses  a  polynomial  basis  for  which  the  generator  [bWb2  -  b2Wb]  is 
not  redundant,  Spencer  (1971).  The  tensor  polynomial  given  above  is  irreducible  and  the  basis  is 
optimal.  Speziale  et  al.  (1991)  have  shown,  for  planar  flows,  that  the  generators  nonlinear  in  the 
anisotropy  tensor  can  be  expressed  in  terms  of  generators  linear  in  the  anisotropy  tensor.  This  fact 
led  to  the  very  simple  form  of  the  SSG  model,  Speziale  et  al.  (1991).  This  fact,  from  a  rigorous 
though  not  necessarily  practical  point  of  view,  limits  the  model  to  planar  flows.  The  results  of 
Speziale  et  al.  (1991)  can  be  used  to  recast  the  present  model  into  its  linear  planar  form.  For 
planar  flows  Speziale  et  al.  (1991)  has  shown  that 

fiipSpj  +  fejpSp,  -  2/3  <  b2S  >  So)  =  -MMw  +  bjrSp,  -  2/3  <  bS  >  6,,)  -  2/3(1  II /b^Sy 

KWpj  +  b%Wpt\  =  -b^[bipWpj  +  bjPWpi} 

from  which  it  follows  that 

KWqpb2}  +  b,pWqpb2qi }  =  (II  +  633633  ){bipWpj  +  bjpWpi] 

and  the  planar  form  of  the  2DMFI  rapid-pressure  model  can  be  written  as 

rtf;  /2<?2  =  (C3  -  2//C3  +  3 1 II<%'  ~  2/3(////633)C7)(5ij] 

+  (C4  -  b33C7){bipSPJ  +  b]VSpi  -  2/3  <bS>  6^}  +  C’"'  <  65  >  [6?  +  211/36,,} 

+  (C5  -  633C8  +  (II  +  b33b33)C9)[b,pWp,  +  b^W^}  +  (C6  <  65  >  +C"'  <  625  >)[6V] 
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which  has  the  same  linear  tensor  bases  as  the  quasi-linear  SSG  model.  The  form  of  the  model 
in  planar  flows  is  consistent  with  to  the  SSG  model  which  is  the  topologically  generic  form  of  a 
general  class  of  models  for  the  fourth  order  rapid- pressure  tensor  Xi}ki ■  The  coefficients  in  the 
model,  however,  are  nonlinear  functions  of  the  invariants  while  in  the  SSG  model,  given  below  for 
comparison, 

IK;  =  -(2C,e  +  C;V)[btJ]  +  C2e[bl  +  2/7/36,,]  +  (C3  -  (-2 11)^)^} 

+  C^k[bxpSP]  +  bjpSpi  -  2/3  <  bS  >  6,,]  +  C$k[bxpWP]  +  bjpWp,\ 

the  C{  are  numerical  constants  and  one  coefficient  in  the  [5,,]  term,  in  the  rapid-pressure  portion  of 
the  model  is  nonlinear.  The  C{  in  the  SSG  model  are  determined  by  a  numerical  optimization  so  that 
the  model  reproduces  as  closely  as  possible:  1)  the  stationary  state  of  the  homogeneous  shear  and 
2)  maximizes  the  kinetic  energy  growth  rate  of  the  rotating  homogeneous  shear  as  close  as  possible 
to  the  Sl/S  =  0.25  predicted  by  rapid  distortion  theory  without  introducing  a  Richardson  number 
similarity,  Speziale  and  Mhuiris  (1989),  while  insuring  that  the  points  of  exchange  of  stability  are 
outside  of  those  predicted  by  the  linear  theory  of  Bertoglio  (1982). 

In  the  return  term  of  the  SSG  model  the  correction  to  the  linear  term  arrived  at  more  or 
less  intuitively  by  Speziale  et  al.  (1991)  is,  to  lowest  order,  vindicated  by  the  present  results. 
Speziale  et  al.  (1991)  have  altered  the  term  linear  in  6,, ,  usually  associated  with  the  slow- pressure 
correlation,  to  include  a  term  involving  the  mean  flow,  a  term  proportional  to  the  production 
V  =  —  2  <  the  usual  — ►  (Ci£  -f  C{V)bi j.  The  present  analysis  indicates  that  the 

portion  of  the  rapid-pressure  contribution  to  the  pressure-strain  correlation,  linear  in  6,j,  has  the 
form  [I?6  <  bS  >  +B4'  <  b2S  >]fv,.  Note  that  <  bS  >  can  be  written  in  terms  of  the  production 
as  <  65  >=  —V/q2.  The  present  analysis  suggests  the  possibility  of  adjusting  the  nonlinear  return 
term  for  mean  velocity  gradient  effects  in  a  similar  way.  In  the  present  model  the  portion  of  the 
rapid-pressure  quadratic  in  the  anisotropy  tensor  is  B+  <  bS  >  [6?  +  2/7/36, ,]  and  therefore  also 
scales  with  the  production. 

Speziale  et  al.  (1992)  have  reflected  on  the  ambiguity  of  the  distinction  between  the  rapid-  and  slow- 
pressure  contributions  to  the  pressure-strain  correlation.  They  have  modeled  the  whole  pressure- 
strain  correlation  but  whether  to  interpret  their  adjustments  (as  described  in  the  previous  para¬ 
graph)  as  incorporating  the  effects  of  the  mean  strain  on  the  return  terms  or  as  contributions  of  the 
rapid-pressure  to  the  total  pressure-strain  is  not  clear.  Here,  however,  the  distinction  is  clear.  The 
present  analysis  starts  from  the  tensor  polynomial  representation  for  the  rapid-pressure  integral 
and  produces  terms  whose  tensor  structure  is  identical  to  that  which  are  traditionally  called  the 
slow-pressure.  The  other  non-linear  pressure-strain  models,  FLT  or  the  SL  model,  Shih  and  Lumley 
(1985),  do  not  have  terms  in  bX}  or  b2}  that  can  be  similarly  identified. 

The  derivation  of  a  rapid-pressure  correlation  satisfying  all  the  limit  states  and  valid  away  from  the 
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limit  state,  though  mathematically  rigorous,  will  not  necessarily  produce  a  model  that  performs 
better  than  other  models  in  flows  for  which  latter  models  are  calibrated.  However,  for  the  general 
class  of  flows  for  which  second-order  models  are  suitable  we  are  of  the  opinion  that  the  satisfaction 
of  all  mathematical  constraints  constitutes  a  necessary  (though  not  sufficient)  requirement  to  assure 
the  predictive  capabilities  in  flows  different  from  those  for  which  the  models  are  calibrated.  The 
existence  of  several  free  parameters  can  then  be  used  to  calibrate  the  model  to  specific  classes  of 
flows  of  computational  interest.  This  tack  is  taken  in  the  next  section  where  the  fixed  points  of  the 
modeled  equation  are  matched  to  the  fixed  points  of  the  homogeneous  shear.  In  this  way  a  model 
for  the  class  of  flows  in  which  the  mean  shear  is  the  dominant  production  mechanism  is  created. 

5.  Calibrating  the  rapid-pressure  correlation  representation 

The  coefficients,  appearing  in  the  2DMFI  model  come  from  first  principles:  they  do  not  result 
from  any  calibration  with  experimental  or  numerical  data.  They  represent  the  minimum  number 
of  determined  coefficients  necessary  to  satisfy  the  mathematical  constraints  on  the  rapid-pressure 
correlation  for  a  three-dimensional  Reynolds  stress.  It,  however,  is  not  a  unique  representation: 
there  are  an  infinity  of  solutions  corresponding  to  different  values  of  the  free  parameters  which,  in 
the  basic  model  representation  shown,  have  been  set  to  zero.  Computations  have  shown  that  the 
predictive  capabilities  of  the  basic  model  are  inadequate.  The  mathematics  built  into  the  model  do 
not  capture  the  experimentally  known  stationary  points  of  homogeneous  shear.  To  compensate  for 
some  of  the  approximations  made  in  the  mathematical  development  the  model  is  modified  so  that 
the  modeled  evolution  equations  have  the  same  asymptotic  behavior  as  that  observed:  to  the  basic 
model  additional  terms  are  added  to  insure  that  the  model  is  consistent  with  an  equilibrium  state 
of  a  particular  benchmark  turbulent  field.  This  is  done  without  sacrificing  any  of  the  mathematical 
principles  built  into  the  model.  It,  however,  can  not  be  done  arbitrarily. 

The  strategy  is  to  require  asymptotic  consistency  with  an  equilibrium  state.  The  modeled  evolu¬ 
tion  equations  are  required  to  have  the  same  fixed  points  as  those  observed  in  experiment.  The 
equilibrium  constraint  is  used  to  obtain  additional  constraints  equations  to  specify  the  free  param¬ 
eters.  This  is  very  similar  to  the  strategy  employed  so  far  in  that  limit  states  are  used  to  set  model 
coefficients.  Here,  of  course,  the  equilibrium  state  is  much  closer  to  those  expected  to  be  seen  in 
flows  of  engineering  interest.  The  free  parameters  will  now  be  called  calibration  coefficients,  A°, 
and  will  be  collectively  denoted  by  X^kl  appearing  in  the  decomposition 

Xuu  =  *§„  +  rx,Fjkl  +  fx°?u. 

Xf?kl  represents  the  additional  constant  terms  necessary  to  capture  the  stationary  state.  This  form 
is  equivalent  to  assuming  that  the  coefficients  in  the  tensor  polynomial  have  the  form  Ai  =  A®  + 
FA?  +  F AI-  The  A\  satisfy  normalization  and  continuity  constraints  and  the  algebraic  equations 
describing  the  stationary  state  of  a  homogeneous  shear.  It  is  at  this  point  that  numerical  or 
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experimental  data  for  the  asymptotic  values  of  the  anisotropy  tensor,  the  production  to  dissipation, 
(P/e) oo,  and  the  ratio  of  time  scales,  (Sk/e) «*,,  are  necessary.  The  details  of  this  strategy  are  now 
given. 

The  A f  satisfy  the  six  equations  given  by  the  homogeneous  form  of  normalization  and  continuity 
constraints,  X^kk  =  0,  X^k  =  0,  which  are  valid  for  all  states  of  the  turbulence.  Substituting  in 
the  calibration  coefficients  allows  six  /If,  to  be  expressed  in  terms  of  seven  free  parameters, 

A\  =  -(50 A%  +  12A§  +  4i4j0)///15  +  (Af,  +  Acu  -  6A13)////5 
A\  =  (20A|  +  3 A|  +  Af0)///15  -  (3 Acn  +  3Af2  +  2A13)////10 
Ac3=  -U/3Al  +  (Acn+3Acl2  +  8A\3)II/3 
A%=  -4/3A%  +  2(A\X  +  2Ai3)///3 
K  =  — 1/3(11A§  +  3A§  +  Af0) 

Ac7=  -2/3(2i4g  +  A\0). 

The  seven  parameters  will  be  determined  by  asymptotic  consistency  with  a  particular  equilibrium 
state.  At  this  point  the  rapid-pressure  correlation  model  is  fully  general  and  it  is  possible  to  write 
it  in  its  final  form,  without  specifying  the  calibration  coefficients,  as 

IK;  /2?2  =  (c3  -  21 1 C3  +  31  lie WSij 

+  e4[bipSPj  +  b]PSp,  -  2/3  <bS>  6ij]  +  C"'  <  bS  >  [6^  +  2///3<5u] 

+  CslbipWp,  +  bJPWvt]  +  [C6  <bS>  +CJ"  <  625  >]6,j 
+  C7[b2pSPJ  +  b%Spt  -  2/3  <  625  >  6tJ] 

+  Cs{b2pWP]  +  b)vWpi\  +  C9[bipWqpb2}  +  b]pWqpb2t}. 

The  calibration  coefficients,  Af,  have  been  added  to  the  base  model  coefficients,  Bt,  to  obtain  the 
final  model  coefficients  C,: 

C3  =  B3-  2F(10Ag  +  3 A%  +  A\0)II/b  -  F(Af,  +  A\2  +  UA\3))III/b 
C3  =  fl£  +  F(A§  +  Af0) 

<T  =  B'^-l/3F(A\l  +  A\2  +  2A\3) 

C4=  B4  +  F(-3Al  +  II(Acn  +  A\2  +  4Acr3)) 

C'”=  B'?  +  F(Acn+A\2  +  4A\3) 

C5=  B5  +  F(-7/3Ac5  +  (~ACU  +  3Af2  +  4Af3)///3) 

C6=  B6  +  F(2A§  +  4Af0) 

C7  =  B7  -  3F(A|  +  Ac9  +  Af0) 

Cs  =  B8-  l/3F(7Ag  +  3A^  -  Af0) 

C9=  B9  +  F(Acn-  A\2) 

Note  that  the  traces  can  be  written  in  terms  of  the  production  <  bS  >=  —V/q2  and  <  b2S  >= 
l/2b33V/q2;  in  planar  flow  <  65  >=  -2bnS  and  <  b2S  >=  -bub^P/q2. 

The  calibration  coefficient  Af  involves  adjustments  to  the  generators  [bS  +  Sb  -  2/3  <  65  >  1] 
and  (bW  -  Wb]  which  are  linear  in  the  anisotropy  tensor.  Numerical  experiments  have  indicated 
that  it  is  important  in  establishing  the  levels  of  the  normal  components  of  the  Reynolds  stresses. 
Combinations  of  Af,  Af  and  Af0  affect  the  generators  quadratic  in  the  anisotropy  tensor,  [b2S  4- 
b2S  -  2/3  <  625  >  1]  and  [b2W  —  Wb2]  and  to  a  very  small  degree  the  S  term.  Combinations  of 
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only  /tg  and  .4,0  can  be  used  to  control  the  contribution  to  f]r  proportional  to  b  and  S.  Experience 
with  the  FLT  model,  Fu  et  al.  (1987)  indicates  that  the  cubic  term  [bWb2  -  b2Wb]  involving 
the  rotation  tensor  is  important  in  controlling  the  relative  level  of  the  normal  components  of  the 
Reynolds  stress  in  situations  with  rotation.  This  suggests  that  the  calibration  coefficients  Aj,,AJ2 
and  v4j3  are  important.  The  combination  of  calibration  coefficients  Acu  -  Aj2  controls  tensor 
products  like  [bWb2  -  b2  Wb]  and  combinations  of  terms  Ajj  ,  A|2  and  A  j3  control  the  contributions 
of  the  b2  and  S. 

The  existence  of  the  free  parameters  allows  the  model  to  be  calibrated  to  a  specific  class  of  flows 
of  computational  interest.  For  example,  in  a  buoyant  flow  one  might  evaluate  the  fixed  points 
using  the  experiments  of  the  homogeneous  shear  in  a  constant  mean  temperature  gradient  thus 
producing  a  model  suitable  for  a  class  of  stratified  flows  of  geophysical  interest.  This  calibration 
would,  of  course,  involve  the  model  for  the  rapid-pressure  correlation  appearing  in  the  heat  flux 
equation.  For  many  engineering  problems  the  primary  production  mechanism  is  the  mean  shear 
and  capturing  the  fixed  points  of  the  homogeneous  shear  in  the  modeled  equations  will  make  the 
model  suitable  for  a  wide  class  of  flows.  At  this  point  one  could  also  consider  using  the  exact 
results  of  rapid  distortion  theory  to  obtain  values  for  the  calibration  coefficients.  This,  however,  is 
inconsistent  with  the  equilibrium  hypothesis  underlying  the  local  approximation  to  the  constitutive 
relation  invoked  for  the  rapid-pressure  correlation  integral.  Moreover  the  rapid  distortion  problem 
is  linear  admitting  a  superposability  that  is  not  possible  in  the  context  of  nonlinear  second-order 
closures.  The  homogeneous  shear  structural  equilibrium  will  be  used  to  fix  the  representation  for 
the  rapid-pressure  correlation. 

The  fixed  points  of  the  homogeneous  shear  are  now  built  into  the  model  by  specifying  the  calibration 
coefficients.  The  modeled  evolution  equations  for  the  homogeneous  turbulence  in  a  mean  velocity 
gradient,  are 

D  j Dt  Ujiij  2c,fcp  K.  UpUj  ^  2tjkp  ^  UpU,  11^  <C  UjUp  ^  u^Up  >  f^!,p 

+  IU  -  Cxebt}  +  C2e[bl  +  2II6tJ/3]  -  2/3eStJ 

D/Dt  £  =  -(Cej  <  uiUj  >  uuj  +  ce2£)  e/k 

assuming  local  isotropy  for  the  dissipation.  The  terms  —C\ebij  -I-  C2£[62_j  +  211  /3bi:]  represent  the 
return  to  isotropy  pressure  correlation.  Using  <  U{Uj  >—  q2(bij  +  1/3 f>a)  the  equations  for  the 
anisotropy  are 

D/Dt  bi:  =  -2(ikpbpjSlk  ~  2<}kpbpiftk  -  [b,pUj,p  +bjpUi,p -2/3£,j  <  bS  >] 

-2/3 Sij  +  2b ij  <bS>+  in,  /q2  ~  (C,  -  2)bt)s/q2  +  C2[&?  +  211 /36^/q2. 

The  mean  strain  and  rotation  tensors  are  StJ  =  I/2S*(<5,|6j2  +  ^12^1)  and  WtJ  =  1/2W‘(6h6}2  - 
6,26ji ),  where,  S  =  U\,2  =  S *  =  W‘ .  Setting  the  D/Dt  =  0  produces  three  algebraic  equations. 
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Inserting  the  experimentally  determined  asymptotic  values  of  btJ  and  Sk/e  produces  three  addi¬ 
tional  constraints  leaving  four  of  the  seven  A\  as  free  parameters.  The  data  from  the  experiments 
of  Tavoularis  and  Corrsin  (1981),  Champagne,  Harris  and  Corrsin  (1970),  Tavoularis  and  Karnik 
(1989),  and  the  DNS  of  Rogers  et  al.  (1986)  are  summarized  in  the  Table  1. 


TC 

CHC 

TK.A 

■iff am 

TK.D 

TK.G 

TK.J 

0.197 

0.137 

0.217 

0.257 

0.197 

0.157 

0.157 

0.147 

0.215 

^2 

-0.14 

-0.165 

-0.165 

-0.165 

-0.17 

-0.148 

-0.154 

-0.149 

-0.158 

-0.143 

-0.083 

-0.133 

-0.143 

-0.133 

-0.113 

-0.103 

-0.093 

-0.153 

^33 

-0.053 

-0.053 

-0.083 

-0.113 

-0.063 

-0.043 

-0.053 

-0.053 

-0.062 

(P/e)  00 

1.75 

1.0 

1.38 

1.37 

1.64 

1.33 

1.45 

1.37 

1.80 

(SK/e)  oo 

6.25 

3.03 

4.2 

4.15 

4.82 

4.5 

4.71 

4.60 

5.7 

Table  1.  Data  for  the  homogeneous  shear  flow. 

There  is  considerable  scatter  in  the  data  due  to  the  technique  used  to  generate  the  turbulence, 
individual  wind  tunnels  in  which  the  different  experiments  were  done,  and  experimental  error.  Not 
all  the  data,  as  has  been  qualified  in  the  references  from  which  the  data  is  drawn,  represent  the 
asymptotic  state.  The  following  fixed  point  values  are  taken  to  be  representative:  =  0.203, 

6j2  =  -0.156,  622  =  -0.143,  633  =  —0.06,  {Sk/e) ^  =  5.54,  (P/e)oo  =  1.73.  They  are  obtained  by  a 
simple  average  of  the  data  of  TC,  TK.D  and  the  DNS.  These  three  cases  are  chosen  because  they 
have  the  highest  values  of  of  the  nondimensional  time  (SK/e)^  corresponding  to  the  flows  that  are 
furthest  in  their  development  to  the  asymptotic  state.  The  values  of  the  invariants  corresponding  to 
these  values  of  the  anisotropy  tensor  are:  //«,  =  -0.057,  ///oo  =  0.0032,  Poo  =  0,573.  Substituting 
these  asymptotic  values  into  the  fixed  point  equations,  D/Dt  bij  =  0,  reduces  the  number  of  free 
parameters  from  seven  to  four. 

Ag  =  -0.29  -  0.06(^o- A|) 

Acn  =  -3.6  +  5^0  -  2Ac13  -  12.74|  -  3.8A§ 

Acu=  -24.5 -  44.24^0- 24^3 +  29A|-8A| 

This  set  of  constraint  equations  is  dependent  on  the  model  for  the  return  to  isotropy  pressure.  For 
simplicity  the  the  nonlinear  return  coefficient  has  been  set  to  zero,  Ci  =  0,  and  the  well  accepted 
value,  C\*  =  3.4,  has  been  chosen.  There  are  still  four  free  parameters.  The  following  values  of  the 
four  free  parameters  are  chosen:  4|  =  0.8, 4|  =  -1.0, 4f0  =  0.01,  A\3  =  0.  The  undetermined  free 
parameters  have  been  set  by  matching  the  values  of  the  anisotropy  for  the  log  layer.  The  procedure 
is  outlined  in  more  detail  in  Appendix  1. 

6.  Computations  and  comparisons  in  homogeneous  turbulence 

The  2DMFI  model  falls  into  the  same  class  of  representations  as  the  FLT  and  SL  models:  they  all 
use  nonlinear  terms  and  invoke  some  form  of  realizability  constraint  to  evaluate  the  coefficients.  For 
this  reason  the  2DMFI  model  will  be  compared  primarily  to  the  nonlinear  SL  and  FLT  models.  For 
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completeness  and  because  it  appears  to  be  a  very  successful  model  for  the  planar  flows,  computations 
with  the  quasi-linear  SSG  model  are  also  shown.  The  SSG  model  is  linear  in  the  anisotropy  tensor 
though  nonlinear  in  that  the  scalar  coefficients  are  functions  of  the  invariants  of  the  the  anisotropy 
tensor.  It  should,  however,  be  kept  in  mind  that  the  SSG  model  satisfies  realizability  for  the  kinetic 
energy  and  not  for  the  individual  Reynolds  stresses  and  is  therefore  in  another  class  of  models.  This 
issue  is  more  fully  explored  in  a  subsequent  section.  Results  are  not  compared  to  the  LRR  model 
as,  in  concordance  with  the  observations  of  Speziale  et  al.  (1991),  the  SSG  model  is  viewed  as  an 
updated  optimized  LRR  model. 

In  all  the  calculations  with  the  2DMFI  model  a  simple  linear  Rotta  type  model  for  the  slow  pressure 
correlation  will  be  used.  This  corresponds  to  C2  =  0  in  the  canonical  form  given  above  and  is 
consistent  with  the  present  calibration  to  the  homogeneous  shear.  For  the  linear  return  coefficient 
a  simple  expression,  C\  —  2  —  31//F1/2,  is  used.  This  satisfies  the  isotropic  limit,  C\  —  2.0,  and 
is  consistent  with  the  assumed  value  for  the  asymptotic  homogeneous  shear,  Cf°  =  3.4.  The  form 
chosen  is  consistent  with  a  weak  form  of  realizability  and  the  recent  results  of  Speziale  et  al.  (1993) 
regarding  the  rate  of  disappearance  of  the  rapid-pressure  correlation  relative  to  the  return  pressure 
correlation  as  the  realizability  limit  is  approached. 

The  values  used  for  the  constants  in  the  dissipation  equation  are:  Cc\  =  1.44,  Ce2  =  1-83.  Note 
that  this  corresponds  to  a  single  universal  fixed  point  (V/e)oo  =  1.88  independent  of  rotation. 
The  single  fixed  point  is  a  well-known  deficiency  common  to  all  the  present  forms  of  the  modeled 
dissipation  equations. 

Case  1:  Homogeneous  shear 

The  calibrated  model  is  now  used  to  compute  the  time  evolution  of  the  homogeneous  shear  flow.  The 
mean  strain  and  rotation  tensors  are  S,j  =  l/2S*(6n6j2  +  ^2^1)  and  WX}  =  1/2VF*(6,i<5j2  —  ^*2^1)1 
where,  S  =  U\, 2=  S*  =  W *.  In  Figure  1  the  time  evolution  of  the  turbulence  energy  is  compared 
to  the  LES  of  Bardina  et  al.  (1983),  and  the  three  models  FLT,  SL  and  SSG.  A  similar  monotonic 
behavior  is  found  for  other  statistics,  612,^11,//  in  the  flow  and,  as  they  do  not  constitute  new  or 
different  information,  are  not  shown.  In  general,  starting  from  physically  realistic  initial  conditions, 
the  flow  attains  its  asymptotic  state  rapidly  and  monotonically.  The  asymptotic  states  which  the 
different  models  attain  aTe  given  in  the  accompanying  table.  The  column  labeled  experimental  data 
is  an  average  of  the  three  cases  TC,  DNS  and  TK.D. 
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Equilibrium 

Values 

2DMFI 

Model 

SL 

Model 

FLT 

Model 

SSG 

Model 

Experimental 

Data 

0.209 

0.202 

0.208 

0.219 

0.203 

*12 

-0.155 

-0.080 

-0.146 

-0.164 

-0.156 

-0.148 

-0.195 

-0.144 

-0.146 

-0.143 

-0.061 

0.007 

-0.064 

-0.073 

-0.06 

manwm 

1.88 

3.42 

1.99 

1.88 

1.73 

(SK/e)  oo 

6.08 

21.35 

6.84 

5.76 

5.54 

Table  2.  Comparison  of  the  model  predictions  for  the  equilibrium  values  in  homogeneous  shear  flow 
(V/e  =  1.88)  with  the  experimental  data  given  in  Tavoularis  and  Karnik  (1989). 

Note  that  the  different  (V/e)^  attained  are  functions  of  the  different  Ce\  and  Ce2  used  in  the  models. 
The  present  form  of  the  dissipation  equation  insures  that  the  quantity  (V/e)^  =  1  88  regardless 
of  initial  conditions  for  all  Sk/e.  This  is  a  shortcoming  of  the  modeled  dissipation  equation  and 
shows  up  in  a  larger  bn  than  the  flow  for  which  it  was  calibrated  in  which  (V/e) <»  =  1 .73 

Case  2:  The  equilibrium  wall  layer 

Another  simple  but  important  test  cases  is  whether  the  model  can  capture  the  stationary  state  of 
the  log-layer  in  the  channel  flow.  The  homogeneous  shear  and  the  log-layer  are  similar  in  that  they 
achieve,  to  a  suitable  approximation,  an  equilibrium  state.  Abid  and  Speziale  (1993)  have  discussed 
the  relevance  of  this  test  case  and  noted  the  inability  of  most  rapid-pressure  closures  to  perform 
successfully  in  the  log-layer.  Our  results  are  in  agreement  with  their  contention  that  a  model  which 
is  asymptotically  consistent  with  the  stationary  states  of  the  homogeneous  shear  will  also  do  well 
in  the  log-layer.  The  models  are  compared  to  the  DNS  of  the  channel  flow  of  Kim  (1993)  which 
is  an  update  of  the  simulations  reported  in  Kim  et  al.  (1987).  The  data  presented  represent  an 
average  of  the  values  of  the  anisotropy  in  the  region  70  <  y+  <  100  outside  the  viscous  sublayer. 


Equilibrium 

Values 

2DMFI 

Model 

SL 

Model 

FLT 

Model 

SSG 

Model 

DNS 

Data 

Experimental 

Data 

bfi 

0.180 

0.079 

0.141 

0.201 

0.180 

0.22 

-0.141 

-0.116 

-0.162 

-0.160 

-0.134 

-0.16 

h°° 

°22 

-0.142 

-0.082 

-0.099 

-0.127 

-0.140 

-0.143 

-0.039 

0.003 

-0.042 

-0.074 

-0.040 

-0.06 

(P/e)  oo 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

(SK/e)  oo 

3.55 

4.30 

3.09 

3.12 

3.73 

3.1 

Table  3.  Comparison  of  the  model  predictions  for  the  equilibrium  values  in  the  log-layer  of  turbulent 
channel  flow  (V/e  =  1)  with  the  DNS  data  of  Kim  (1993)  and  the  data  of  Laufer  (1951)  given  in 
Abid  and  Speziale  (1993). 


Case  3:  Homogeneous  shear  with  rotation 

The  present  test  case,  the  homogeneous  shear  with  rotation,  and  the  next  test  case,  the  homogeneous 


shear  with  streamline  curvature,  are  important  test  cases.  In  both  these  cases  additional  forces, 
which  stabilize  or  destabilize  the  flow,  are  present.  These  effects  appear  in  the  evolution  equations  as 
additional  production  mechanisms  for  the  Reynolds  stresses.  In  the  case  of  rotation  the  production 
terms  in  the  evolution  equation  for  the  turbulence  kinetic  energy  do  not  directly  depend  on  the 
rotation:  the  turbulence  energy  production  depends  on  the  rotation  only  through  the  off-diagonal 
components  of  the  Reynolds  stress  while  in  the  case  of  streamline  curvature  the  production  terms 
in  the  evolution  equation  of  the  turbulence  kinetic  energy  do  directly  dependent  on  the  curvature. 
These  two  cases  are  important  test  cases  not  only  because  1 )  the  models  have  not  been  calibrated  for 
them  but  also  because  2)  the  model  will  have  to  predict  both  the  stabilization  and  the  destabilization 
of  the  turbulence  and  3)  the  critical  values  of  the  governing  parameters  which  demarcate  the  regions 
of  flow  stabilization  from  flow  destabilization. 

For  flows  in  the  rotating  coordinate  system  the  Coriolis  terms  must  be  carried  and  Wv  appearing 
in  Ho  must  be  replaced  by  the  total  rotation  tensor  VFtJ  -Ffjifcftfc-  Thus  in  the  rapid-pressure  model 
W  =  S(l-2Q/S). 

Figure  2b-2e  show  how  the  models  perform  in  rotating  shear  for  rotation  to  shear  ratios  n/S  = 
0,  1/4,  1/2  compared  to  the  LES  data  of  Bardina  et  al.  (1983)  in  Figure  2a.  In  general,  all  the 
models  are  able  to  capture  both  flow  stabilization  for  some  A  >  Q/S  >  B  and  flow  destabilization 
for  some  A  <  Q/S  <  B.  The  points  A  and  B  represent  the  points  of  neutral  stability  on  a 
bifurcation  diagram  in  the  phase  plane  ( e/Sk )<*>  and  ( ft/S )(X>.  All  of  the  models  have  a  bifurcation 
diagram  of  the  same  general  form,  Speziale  and  Mhuiris  (1989),  indicating  a  stabilization  of  the 
flow  outside  of  some  region  of  approximate  size  0  <  fi/5  <  0.5,  predicted  by  the  linear  rapid 
distortion  theory  of  Bertoglio  (1982).  The  most  important  facts  concerning  the  different  models 
for  the  homogeneous  rotating  shear  can  be  summarized  by  indicating  the  unstable  regions  in  which 
the  models  predict  a  non-trivial  equilibrium  (e/Sk) 


SSG : 

-0.09  < 

n/s  <  0.53 

RDT  : 

0.00  < 

n/s  <  o.5o 

2 DM  FI  : 

-0.063  < 

n/S  <  0.502 

SL  : 

-0.14  < 

n/s  <  o.4o 

FLT  : 

-0.11  < 

n/S  <  0.39 

Near  the  point  of  linear  neutral  stability  f l/S  =  1/2  both  the  SL  and  the  FLT  models  predict  a 
premature  restabilization  at  values  of  fi/5  20%  and  22%  lower  than  that  predicted  by  the  linear 
theory.  The  2DMFI  model  is  within  4%  of  the  linear  prediction. 

None  of  the  models  tested,  linear  or  nonlinear,  captures  the  point  of  maximum  kinetic  energy 
growth  at  (Q/ S)mriT  =  0.25.  To  do  so  would  mean  that  the  equations  would  exhibit  a  Richardson 
number  similarity  which,  as  Speziale  and  Mhuiris  (1989)  have  shown,  is  not  admitted  by  the  Navier 
Stokes  equations.  The  two  models  that  come  closest  to  (Cl/S)max  =  0.25  are  SSC1  at  fi/5  =  0.22. 


20 


which  was  calibrated  using  this  fact,  and  2DMFI  at  f l/S  =  0.2  which  was  not  calibrated  using 
any  rotating  flows.  The  current  modeled  dissipation  equations  predict  a  ( V/e )oo  =  -26i25Jfc/e  = 
(Ce 2  —  l)/(Cei  —  1)  =  const.,  where  the  constant  is  model  dependent  but  independent  of  rotation 
rate.  The  constant  attains  the  value  for  the  equilibrium  homogeneous  shear  for  arbitrary  rot.  tion 
rate,  a  fact  which  is  not  consistent  with  observation.  In  a  flow  that  is  stabilized  by  rotation,  say 
ft/5  =  1,  production  must  be  less  than  dissipation  for  the  equilibrium  state  to  be  reached.  The 
dissipation  equation  cannot  be  used  for  calibration  in  rotating  flows  without  compromising  the 
model  when  a  dissipation  equation  capable  of  predicting  the  stationary  values  of  (V/e)^  =  /(ft/5) 
becomes  available. 

Case  4:  Homogeneous  shear  with  curvature 

For  the  homogeneous  shear  with  streamline  curvature  the  mean  strain  and  rotation  tensors  are  StJ  = 
l/25*(Mj2  +  M,'i)  and  Wi}  =  \/2W’{StXb]2-6i26]X)  where  S’  =  S(l-stb)  and  W’  =  5(1  +stb) 
where  stb  =  ( Uc / Rc)/S  is  the  stability  parameter.  The  geometry  for  the  curved  homogeneous  shear 
follows  that  of  Holloway  and  Tavoularis  (1992):  Rc  is  the  radius  of  curvature  of  the  flow  and  Uc 
is  the  axial  velocity  at  the  centerline  and  the  crosstream  gradient  of  the  axial  velocity  is  the  shear 
Ui,2  =  5.  The  stability  parameter  has  been  renamed  stb  so  as  not  to  confuse  it  with  5  which  is 
traditionally  used  for  the  mean  shear.  The  kinetic  energy  growth  rate  is  suppressed,  relative  to  the 
homogeneous  shear,  for  stb  >  0  and  increased  for  stb  <  0,  while  for  stb  >  0.05  the  experimental 
data  indicates  a  relaminarization. 

Figure  3  compares  the  model  results  to  the  experimental  data  for  6J2  versus  stb.  The  plot  has  been 
generated  by  computing  the  flow  from  the  beginning  of  the  straight  section  of  the  wind  tunnel  to  a 
value  of  St  =  10  which  corresponds  to  the  end  of  the  curved  section.  The  initial  conditions  on  the 
second  order  moments  are  given  by  the  experimental  data.  The  initial  condition  on  the  dissipation 
rate  is  determined  by  matching  to  the  kinetic  energy  growth  rate  at  the  beginning  of  the  straight 
section. 

The  different  models  all  capture  the  trend  in  the  stabilization/destabilization  with  respect  to  the 
stability  parameter.  The  primary  difference  in  the  predictions  of  the  different  models  seems  related 
to  their  ability  to  capture  the  homogeneous  shear  at  stb  =  0.  The  results  of  the  different  models 
would  be  in  more  agreement  for  negative  and  small  positive  values  of  stb  if  they  predicted  the  same 
results  for  the  homogeneous  shear. 

The  bifurcation  diagram  for  the  second  order  models,  in  Figure  4,  was  generated  by  letting  the 
solution  procedure  go  to  its  asymptotic  state.  There  is  a  critical  value  of  stbc  at  which  the  stabi¬ 
lizing  effects  of  curvature  begins  to  causes  a  negative  kinetic  energy  growth  rate,  which  ultimately 
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relaminarizes  the  flow.  The  critical  values,  stbc,  predicted  by  the  different  models,  are 

HkT  :  stbc  =  0.05 
2DMFI  :  stbc  =  0.067 
FIT  :  stbc  =  0.075 
SSG  :  stbc  =  0.10 
SL:  stbc  =  0.105 

where  H kT  is  from  the  experimental  data  of  Holloway  and  Tavoularis  (1991 ).  There  is  a  consistent 
trend  for  the  SSG,  2DMFI  and  FLT  models,  when  compared  to  the  critical  values  for  the  rotating 
shear:  the  higher  (Sl/S)c  for  stabilization  of  the  flow  correspond  to  higher  stbc.  The  SL  model  has 
small  but  nonzero  e/Sk  over  the  range  0.25  <  stb  <  stbc. 

Case  5:  Two  and  three-dimensional  strains 

The  rapid  pressure  model  is  used  to  compute  three  strain  flows:  the  plane  strain  and  the  axisym- 
metric  contraction  and  expansion.  These  flows  are  another  test  case  as  the  rapid  pressure  model 
has  not  used  these  flows  to  set  the  calibration  coefficients.  The  results  are  compared  to  the  direct 
numerical  simulations  of  Lee  and  Reynolds  (1985).  Because  the  simulations  are  conducted  at  low 
Reynolds  number  the  anisotropies  are  expected  to  be  somewhat  higher  than  those  of  a  fully  de¬ 
veloped  turbuelence.  However  the  use  of  the  physical  experiments  conducted  at  higher  Reynolds 
number  is  also  somewhat  tenuous  as  the  initial  conditions  on  e/Sk,  as  has  been  pointed  out  by 
Speziale  et  al.  (1991),  are  not  known  with  certainty.  The  same  test  cases  as  those  given  in  Speziale 
et  al.  (1991)  are  used.  Our  results  are  also  compared  to  the  SSG  model  as  it  appears  to  be  the 
current  model  that  gives  the  best  results.  The  evolution  of  the  kinetic  energy  for  these  flows  is  not 
presented;  the  results  for  the  models  are  in  very  good  agreement  with  the  data  and  each  other  and 
do  not  definitively  distinguish  between  the  various  models. 

Figure  5  shows  the  evolution  of  the  anisotropy  for  the  plane  strain,  Sjj  =  S*($»ifyi  —  fofyj) 
starting  from  isotropic  initial  conditions.  Figure  6  and  7  shows  the  evolution  of  the  anisotropy  for 
the  axisymmetric  contraction  and  expansion.  Here  for  the  contraction  Sij  =  -  l/26<2^2  - 

1/26*3^3).  For  the  expansion  S  is  replacesd  with  -5.  Results  for  all  the  plane  and  axisymmetric 
strain  flows  capture  the  trends  nicely. 

8.  General  discussion 

The  constraints  of  geostrophy,  realizability,  joint-realizability,  normalization  and  continuity  have 
been  used  to  create  the  2DMFI  model.  The  use  of  a  realizability  type  constraint  to  obtain  values 
for  unknown  coefficients  in  the  models,  in  the  case  of  the  strong  form  of  realizability,  has  been 
criticized  on  the  grounds  that  one  should  not  use  an  extreme  state  to  set  the  coefficients  in  an 
equilibrium  model.  This  certainly  is  the  case  in  any  rigorous  interpretation  of  the  statement  - 
requiring  all  the  scales  of  the  motion  to  satisfy  the  indicated  limit.  However,  from  the  point  of 
view  of  a  useful  engineering  approximation,  it  should  be  kept  in  mind  that  <  utu;  >  represents 
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an  integral  over  all  scales  of  the  motion,  from  the  production  scales  at  nt  ~  1  and  larger  to  the 
dissipation  scales  kt}  ~  1.  In  a  turbulence  with  a  k~5/3  inertial  subrange  in  which  there  is  enough 
of  a  separation  of  scales  for  a  second-order  simulation  to  be  useful,  say  at  least  Re/  ~  104,  the  ratio 
between  the  dissipative  and  the  energy  containing  length  scales  is  t]/(  ~  ReJ3^4  ~  1000  and  the 
flow  scales  range  over  0  <  kI  <  1000.  However,  approximately  8551  of  the  energy  of  the  motion  is 
contained  in  the  first  decade  nt  <  10:  the  major  contribution  to  <  u,  u3  >  is  from  the  scales  of  the 
motion  greater  than  one  tenth  of  the  production  scales,  for  which  K/i  ~  1.  If  only  the  largest  1% 
of  the  flow  scales,  i.e.  from  0  <  k(  <  10,  begins  to  lose  an  eigenvalue  of  the  Reynolds  stress  tensor, 
through  some  dynamical  or  kinematical  agency,  <  utiij  >  begins  to  approach  the  realizable  limit. 

Part  of  the  reluctance  to  accept  the  second-order  modeling  technique  is  that  turbulence  simulations 
of  complex  inhomogeneous  flows  with  multi-dimensionai  mean  flows  which  may  have  body  forces, 
streamline  curvature  or  rotational  effects,  are  extremely  difficult  to  compute.  ^From  a  strictly 
practical  point  of  view  -  computability  -  incorporating  the  realizability  constraints  into  the  models 
for  unknown  correlations  has  some  very  tangible  and  beneficial  efferts.  During  the  convergence  to 
a  solution,  from  more  or  less  arbitrary  initial  conditions,  the  iteration  will  be  plagued  with  negative 
normal  stresses  and  correlation  coefficients  larger  than  one.  When  this  occurrs  the  solution  is  clipped 
and  the  solution  procedure  restarted  from  the  new  clipped  initial  conditions.  The  frequency  of  this 
clipping  and  resetting  procedure  is  substantially  reduced  (in  simple  flows  even  eliminated)  when 
using  realizable  models. 

How  one  gets  to  a  realizable  turbulence  is  not  an  issue  in  problems  with  steady  states,  as  long  as 
the  final  state  is  realizable.  This,  however,  is  not  the  case  for  problems  that  are  unsteady.  For  time- 
varying  flows,  varying  on  the  integral  time-scale,  the  quasi-steady  problem  for  which  the  second- 
order  methods  are  still  suitable  -  this  issue  becomes  serious.  Our  experience  with  the  Reynolds 
averaging  procedure,  in  buoyantly  driven  elliptic  flows  with  rotation,  indicates  that  it  acts  as  a  low- 
pass  filter:  the  rapidly  fluctuating  instantaneous  dynamics  are  subsumed  by  the  averaging  procedure 
leaving  the  slow-time  large-scale  parts  of  the  flow,  evolving  on  time  scales  commmensurate  with  the 
integral  time  scale,  to  be  captured  by  the  computation.  If  the  simulation  is  to  reflect  the  physics 
of  the  time  evolution  of  the  flow,  it  must  stay  realizable.  Clearly  excessive  realizability  violations, 
requiring  a  clipping  and  resetting  of  the  solution,  which  produces  a  solution  that  never  evolves  far 
from  the  transient  associated  with  the  most  recent  clipped  initial  condition,  are  not  acceptable. 
In  such  flows,  satisfying  the  realizability  constraint  has  very  important  consequences  regarding  the 
validity  of  the  time  evolution  of  the  flow. 

The  coefficients  of  the  tensor  polynomials  used  to  represent  an  unknown  correlation  in  a  constitutive 
relation  are,  according  to  the  theory,  non-constant  functions  of  the  invariants  of  the  independent 
tensors  and  thus  depend  on  the  state  of  the  turbulence.  These  coefficients,  in  "realizable”  turbulence 
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models  using  some  form  of  the  strong  realizability  constraint,  are  obtained  by  requiring  that  the 
rate  of  change  of  a  positive  semi-definite  quantity  be  zero  when  some  limit  state  is  achieved.  The 
coefficients  so  obtained  are  constants  and  are  strictly  valid  only  at  the  limit  state.  To  the  limit  state 
value  of  the  coefficients,  which  are  regarded  as  a  zeroeth  order  approximation  to  the  coefficient  in 
a  general  flow,  ad  hoc  corrections  are  then  added.  These  corrections,  which  vanish  as  the  limit 
state  is  approached,  require  some  sort  of  numerical  optimization  with  well  known  simple  flows. 
This  summarizes  the  methodology  used  in  other  "realizable”  turbulence  models.  In  the  linear  SSG 
model  the  coefficients  are,  for  the  most  part,  constants  that  are  also  set  by  matching  to  a  limit  state. 
In  the  case  of  the  SSG  model  the  structural  equilibrium  limit  is  used  to  set  the  constants  in  the 
model.  This  is  a  much  closer  approximation  to  the  turbulence  expected  to  be  seen  in  engineering 
problems.  The  constants  in  the  SSG  model  may  be  viewed  as  the  values  the  nonconstant  coefficients 
have  near  the  equilibrium  state. 

In  the  present  method  the  realizability  principles  are  used  to  obtain  non-constant  values  of  the 
coefficients  in  the  constitutive  relations  that  are  valid  for  all  states  of  the  mechanical  turbulence. 
Recall  that  the  basic  form  of  the  rapid-pressure  model  is  comprised  of  two  parts  Xfjkl  and  FX?kl, 
the  ’’off- realizability”  correction.  X®kt  satisfies  simultaneously  the  five  constraints  -  geostrophy, 
realizability,  joint-realizability,  normalization  and  continuity,  while  FX[]kl ,  also  obtained  analyt¬ 
ically,  satisfies  the  substantially  less  extreme  joint-realizability  as  well  as  the  homogeneous  form 
of  the  normalization  and  continuity  constraints.  Thus,  although  the  model  is  consistent  with  the 
realizable  limit  state,  the  values  of  the  model  coefficients  are  not  the  values  the  coefficients  have  at 
the  limit  state.  An  extreme  state  of  the  flow  has  only  been  used  to  set  the  coefficients  in  the  Xfjkl 
part  of  the  model  -  the  additional  FXf]kl  terms  are,  for  the  mechanical  turbulence,  fully  general. 

It  is  to  this  basic  model,  valid  for  all  states  of  turbulence  not  just  at  the  realizability  limit,  that 
one  adds  the  FX^k,  term  that  is  necessary  to  obtain  the  equilibrium  values  in  the  flows  evolution 
from  arbitrary  initial  conditions.  The  requirement  of  asymptotic  consistency  with  an  equilibrium 
state,  first  used  by  Speziale  ct  al.  (1991),  is  the  single  most  important  and  consistent  physical 
requirement  one  can  impose.  Second-order  closure  methodology  is  built  around  the  assumption 
that  there  is,  to  a  suitable  approximation,  for  the  class  of  flows  to  which  second-order  methods 
are  appropriate,  a  equilibrium  state  and  in  the  absence  of  disturbing  forces  the  flow  relaxes  to 
that  state  on  a  time  scale  similar  to  the  eddy  turnover  time.  It  is  precisely  this  phenomenological 
behavior  that  is  built  into  the  model  by  requiring  the  fixed  points  of  the  modeled  equations  to  be 
consistent  with  those  obtained  from  experiment.  The  assumption  that  allows  the  parameterization 
of  the  two-point  correlation  as  a  local  function  of  the  anisotropy  tensor,  Xiju  =  Xijki(b{j),  is  such 
an  equilibrium  assumption. 

The  penalty  paid  for  these  additional  features  associated  with  the  satisfaction  of  the  mathematical 
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constraints  is  a  more  complex  model.  It  should,  however,  be  pointed  out  that  the  present  2DMFI 
model  has  the  same  tensor  basis  as  the  FLT  model  and  is  therefore  no  more  complex  except  for 
the  expressions  for  the  nonconstant  coefficients.  Moreover,  the  penalty  is  slight  in  the  light  of 
the  reduction  of  the  computational  difficulties  found  during  the  calculation  of  quasi-steady  time 
evolving  flows  with  this  representation  for  the  rapid-pressure.  The  model,  along  with  several  other 
models,  has  been  used  to  compute  inhomogeneous  buoyancy  driven  rotating  flows  that  occur  in 
the  Czochralski  crystal  growth  melt  in  which  the  Reynolds  stress  are  three-dimensional,  Ristorcelli 
and  Lumley  (1991a,  1993),  Ristorcelli  (1991).  In  computing  these  time  varying  flows  it  was  found 
that  the  present  2DMFI  model  produced  far  fewer  realizability  violations  during  the  course  of  the 
flows  evolution.  For  a  quasi-steady  flow  this  is  a  crucial  point:  every  time  realizability  is  violated 
the  solution  is  reset  and  the  solution  never  evolves  past  the  transients  associated  with  reseting  the 
initial  conditions.  Such  a  computation  cannot  be  expected  to  reflect  an  ensemble  average  of  the 
original  system. 

9.  Suggestions  for  future  work 

In  the  effort  to  produce  a  representation  for  the  rapid- pressure  correlation  valid  for  the  class  of 
flows  to  which  second-order  modeling  is  suitable,  some  shortcomings  in  the  data  on  homogeneous 
"building  block”  flows  have  become  apparent.  Though  the  homogeneous  shear  seems  reasonably 
well-documented  it  is  not  clear  that  the  asymptotic  states  have  been  reached  in  some  of  the  exper¬ 
iments.  Moreover  the  discrepancy  between  the  high  values  of  obtained  in  the  DNS  versus  those 
seen  in  the  laboratory  data  has  not  been  explained.  Additional  work  expanding  on  the  notion  of 
two  classes  of  flows,  according  to  Tavoularis  and  Karnik  (1989),  and  therefore  the  possibility  of  two 
equilibrium  states,  merits  investigation. 

For  the  homogeneous  shear  with  rotation,  a  very  basic  flow,  there  doesn’t  seem  to  be  substantial 
data,  LES,  DNS  or  experimental  definitively  describing  its  stationary  states.  At  the  very  least  an 
assessment  of  the  bifurcation  diagram  (e/Sk) «,  versus  (0,/S)oo  predicted  by  the  linear  theory  would 
be  useful.  Also  the  equilibrium  values  (V/e)oo  would  be  useful  for  further  developments  regarding 
the  dissipation  equation’s  dependence  on  rotation.  This  may  remedy  the  under-prediction  of  the 
kinetic  energy  growth  rates  as  a  function  of  Sl/S  for  all  the  models.  The  present  class  of  dissipation 
equations  predicts  an  asymptotic  state  in  which  (P/e)oo  is  a  model  dependent  constant,  which  for 
all  rotation  rates  has  the  same  value  that  it  does  in  the  asymptotic  shear.  Had  the  stationary 
values  of  the  anisotropy  tensor  and  (V/e)oo  and  (e/Sk)0 0  been  available  application  of  the  present 
methodology  would  have  produced  a  set  of  modeled  evolution  equations  whose  fixed  points  matched 
to  the  fixed  points  of  the  rotating  shear,  independent  of  the  deficiencies  in  the  dissipation  equation. 

The  present  rapid-pressure  model  is  expected  to  distinguish  itself  in  complex  three-dimensional 
flows.  For  the  planar  flows  for  which  test  cases  exist  the  model  out-performs  the  nonlinear  models 
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using  realizability  type  constraints.  It  is  only  moderately  better  than  the  topologically  generic 
form  of  the  SSG  model  suitable  for  simple  planar  flows.  It  is  unfortunate  that  there  does  not  exist 
any  suitable  DNS  or  LES  of  flows  in  which  the  presence  of  a  body  force  causes  the  larger  scales 
of  the  motion  to  become  quasi-two-dimensional.  Such  a  test  case  would  help  further  establish  the 
utility  of  incorporating  some  of  the  more  complex  physics  into  the  structure  of  the  model  as  well 
as  perhaps  pointing  out  the  potential  deficiencies  of  models  developed  using  two-dimensional  mean 
flows. 

10.  Summary  and  Conclusions 

A  new  representation  for  the  rapid-pressure  strain  correlation  with  a  minimum  of  ad  hoc  constants 
has  been  devised.  The  rapid-pressure  model  produces  the  proper  behavior  in  five  different  limits: 

1)  the  geostrophic  limit  in  which  the  eigenvalue  of  the  Reynolds  stress  tensor,  <  >,  along  the 

axis  of  rotation  vanishes, 

2)  the  realizable  limit  in  which  an  arbitrary  eigenvalue  of  the  Reynolds  stress  vanishes, 

3)  the  joint-realizable  limit  in  which  an  eigenvalue  of  <  Od  ><  u,Uj  >  -  <  6u,  ><  du,  >  vanishes, 

4)  the  isotropic  limit  in  which  the  anisotropy  tensor,  6,y  =  0, 

5)  the  asymptotic  limit  of  a  structural  equilibrium  in  which  D/ Dt  bt]  =  0. 

The  model  has  the  general  form 

Xm  =  x°u  +  Fxru  +  fx*u 

where  the  Xijki  are  polynomials  in  the  anisotropy  tensor.  X°kl  satisfies  the  five  constraints:  the 
limit  states  of  geostrophy,  realizability,  joint-realizability,  and  the  integral  constraints  of  continuity 
and  normalization.  X[jkl  satisfies  the  three  constraints  of  joint-realizability,  continuity  and  normal¬ 
ization.  Both  Xfjkl  and  FX[jkl  are  obtained  analytically:  they  represent  the  simplest  analytical 
expressions  that  are  capable  of  satisfying  all  the  mathematical  constraints.  X°°kl  is  obtained  by 
requiring  asymptotic  consistency  with  the  structural  equilibrium  state  for  homogeneous  shear.  Ex¬ 
perimental  data  is  required  to  determine  X°°kl .  The  stationary  values  of  the  anisotropy  tensor,  b°°, 
and  (T’/s)00  and  (s/Sk)x  are  inserted  into  the  modeled  evolution  equations  to  obtain  the  structural 
equilibrium  component  of  the  rapid-pressure,  This  insures  that  the  fixed  points  of  the  mod¬ 

eled  equations  match  the  experimentally  determined  fixed  points.  This  process  has  been  carried 
out  for  the  homogeneous  shear  producing  a  model  that  should  be  accurate  for  three-dimensional 
mean  flows  in  which  the  mean  shear  is  the  predominant  production  mechanism. 

Several  points  regarding  the  2DMFI  model  merit  mention: 

1)  All  the  coefficients,  A,,  in  the  X°kl  +  FXfjkl  portion  of  the  model  are  obtained  from  first 
principles  and  are  valid  for  all  states  of  the  mechanical  turbulence.  They  take  on  their  values  at 
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the  limit  states  only  when  the  limit  states  are  achieved  unlike  other  models. 

2)  The  present  form  of  the  model  is  consistent  with  the  stationary  state  of  homogeneous  shear.  This 
insures  that  in  the  absence  of  any  disturbing  forces,  the  predicted  flow  will  relax  to  the  equilibrium 
state. 

3)  The  model  is  at  present  the  only  model  which  is  consistent  with  predictions  of  the  Taylor- 
Proudman  theorem  assuring  that  the  modeled  equations  are  frame-indifferent  when  the  components 
of  the  Reynolds  stress  along  the  axis  of  rotation  vanish.  This  frame  invariance  is  the  most  notable 
feature  of  the  model  and  is  expected  to  be  important  for  the  computation  of  engineering  and 
geophysical  flows  in  which  body  forces  play  an  important  role.  In  meteorological  flows  the  present 
2DMFI  rapid -pressure  correlation  is  consistent  with  the  geostrophic  limit  atmosphere  attains  above 
the  planetary  boundary  layer.  As  such  the  model  will  be  relevant  to  the  computation  of  mesoscale 
meteorological  flows. 

It  is  expected  that  the  satisfaction  of  the  2DMFI  principle  will  be  important  in  three-dimensional 
flows  in  which  the  largest  scales  of  the  motion  are  two-dimensionalized  by  body  forces  or  kinematic 
constraints.  These  flows  include:  1)  turbulence  in  which  a  strong  stable  stratification  suppresses 
the  vertical  component  of  the  velocity  field;  2)  turbulence  in  the  flows  affected  by  magnetic  fields 
3)  turbulence  influenced  by  rotation  such  as  those  in  crystal  growth  processes  or  occurring  in  tor¬ 
nadoes,  swirl  combustois  and  turbines;  4)  turbulence  influenced  by  Coriolis  forces  in  which  the 
largest  scales  of  the  flow  undergo  a  Taylor- Proudman  reorganization  ( Ro  <  2)  as  might  occur  in 
large  scale  geophysical  flows;  5)  turbulence  near  a  free  surface  at  which  one  of  the  components  of 
the  fluctuating  velocity  is  suppressed  and  the  mean  flow  is  two-dimensional;  6)  the  environmen¬ 
tally  important  shallow  water  free  shear  flows  as  the  shallow  water  jet  associated  with  waste  heat 
exchangers,  near-shore  pollution  dispersal,  tidal  estuary  flows,  and  mixing  associated  with  thermal 
and  salinity  inflows  in  lakes  and  rivers.  However,  until  suitable  data  bases,  DNS  or  LES,  of  these 
complex  flows  with  body  forces  become  available  the  full  potential  of  a  rapid-pressure  model  built 
from  first  principles  in  three-dimensional  flows  can  not  be  verified.  The  present  model  does  however 
reproduce  the  experimental  data  at  least  as  well  as  the  currently  available  models  for  a  wide  class 
of  two  dimensional  flows. 
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Appendix  1:  A  synopsis  of  the  final  rapid-pressure  correlation  representation 

For  convenience  and  clarity  the  final  form  of  the  model  is  summarized  here.  The  general  final  form 
of  the  model  is: 


Ho  /2 q2  =  [C3  -  2//C3  +  SIIIC^S,, 

+  C4[biPSPJ  +  bjpSpi  -  2/3  <bS  >  +  C'"  <  bS  >  [6?  +  277/3.5,,] 

+  CsKWpj  +  bJPWpi]  +  [C6  <bS>  +C”  <  b2S  >}bv 
+  CrKSpj  +  b^Sp,  -  2/3  <  b2S  >  Sv] 

+  C8[b?pWPJ  +  bjpWpt]  +  C9[b,pWqpb2q3  +  b]pWqpb2qt}. 

For  flows  in  rotating  coordinate  systems,  WX]  appearing  in  fll,  must  be  replaced  by  the  total  rotation 
tensor  W{j  +  The  coefficients  C,  are  given  by  a  sum  of  the  basic  model  coefficients,  Bt, 

which  come  from  first  principles,  and  the  calibration  coefficients  A'.  The  basic  model  coefficients 
are  related  to  the  coefficients  A,  in  the  fourth  order  tensor  polynomial  expression  by 

B3  =  2(A,  +  A2)  =  (2/27)[41  +  42/7  -  0.1F(221  +  42077)]/77d 
B'{  =  A9  +  >l,o  =  -( 14/3)(  1  -I-  377)/77d  +  0.6F/(  1  +  377) 

B'f  =  -1/3(A„  +  A,2  +  2,4,3)  =  (55  +  8477)/377<, 

B4  =  A3  +  Aa  +  2 As  =  3/7/d  -  0.9F/(1  +  377) 

B'H  =  An  +  A,2  +  4A,3  =  -9/ 1  Id 
Bs=  A3~A4  =  — ( l/30)(  10  +  21F)/(1  +  377) 

B6  =  2  A9  +  4A,0  =  -18II/IId  +  3F/(1  +  377) 

B7  =  A6  +  A7  +  2A8  -  2A9  -  2A,0  =  ~9lIId  -  1.8F/(1  +  377) 

Bg  =  A6~  A7  =  (i/5)(3F—  5)/(l  +  377) 

A„ -A, 2  = -3/(1  +  377) 

The  A,  are  given  in  Appendix  3.  The  Ci,  which  reflects  the  application  of  the  additional  constraint, 
asymptotic  consistency  with  an  equilibrium  state,  are  obtained  from  the  77,  by  C,  =  Bi  +  F Ac%. 


They  are  given  by 

c3  = 

B3  -  2F(10A|  +  3 Ac9  +  A,0)/7/5  -  F(A,,  +  A\2  +  14AJ3))777/5 

B3  +  F(A|  +  A^0) 

c3  = 

B3  -  1/3F(A5,  +  A^j  +  2Af3) 

c4  = 

B4  +  F(-3Ag  +  77(A,X  +  Af2  +  4Af3)) 

II 

17 

7?"'  +  F(Af,  +  Af2  +  4Af3) 

c5  = 

B-,  +  n-7/3A|  +  (-A,,  +  3A^2  +  4A^)77/3) 

c6  = 

fi6  +  F(2A|  +  4Af0) 

c7  = 

Bt  -  3F(A|  +  A$  +  Af0) 

c8  = 

Bg  -  l/3F(7Ag  +  3A9  -  A,0) 

c9  = 

F9  +  F(Af,-Af2) 

where  the  A£  are  expressed  in  terms  of  the  seven  free  parameters  A£{5,8  -  13}.  Without  further 
specifying  the  calibration  coefficients,  A^,  the  above  rapid-pressure  model  is  general,  suitable  for  the 
flows  for  which  second-order  modeling  is  suitable.  Calibration  to  a  particular  archetypal  flow,  by 
matching  the  fixed  points  of  the  modeled  evolution  equations  with  experimentally  determined  fixed 
points,  will  result  in  a  model  suitable  for  diverse  flows  within  that  class  of  flows.  For  a  flow  in  which 


31 


the  dominant  production  mechanism  is  associated  with  the  mean  shear  the  hommogeneous  shear 
is  used  to  set  the  model  coefficients.  Asymptotic  consistency  with  the  homogeneous  equilibrium 
shear  and  a  linear  model  for  the  return  coefficient  produce  the  following  values  for  the  calibration 
coefficients: 

A%  =  -0.29  -  0.06(Aio  -  A§) 

Ajj  =  -3.6  +  5Af0  -  2 Acl3  -  12.7 Ag  -  3.8 A| 

A\ 2  =  -24.5  -  44.2Af0  -  2 A{3  +  29A|  -  8A| 

where  Ag  =  0.8.  A§  =  -1.0,  A\Q  =  0.01,  Aj3  =  0. 

The  homogeneous  form  of  the  continuity  and  normalization  conditions  have  been  used  to  produce 
relations  between  the  free  parameters.  They  are: 

A\  +  4 Ac2  -  2Ac8II  +  III(ACU  +  A\2  +  2Af3)  =  0 
^3  +  ^4  +  5A1  -  II{A\X  +  AX2  +  4A\3)  =  0 
A%  +  Ac7  +  5  A%  +  A§  +  Af0  =  0 

3Af  +  2 Ac2  -  2Ag/7  +  4Aj3//7  =  0 
3A$  +  4 Al  -  277(Af,  +  2Aj3)  =  0 
3A7  +  4  Ag  +  2Aj0  =  0 

The  fact  that  the  A£  must  satisfy  these  equations  has  been  used  to  express  six  of  the  A^{1  -  4,6, 7} 
in  terms  of  the  seven  Af{5,8  -  13}  in  the  expressions  for  the  C,  =  +  FA 

Appendix  2:  The  calibration  coefficients  in  the  rapid-pressure  representation 

The  model  has  been  calibrated  to  match  the  fixed  points  of  the  homogeneous  shear.  The  calibration 
process  is  described  more  detail  in  this  appendix.  The  following  equations  describe  the  evolution 
of  the  anisotropy: 

D/ Dt  btj  =  -2(,kpbP]ilk  -  2(jkpbp,ilk  -  [bipUj,p  +bjpUi,p  -2/3 6^  <  bS  >] 

-2/35, j  +  26, 2  <bS>+  n:,  /q2  ~  (Ci  -  2)bt]e/q2  +  C2[b?}  +  2II/36ij}e/q2 

For  the  planar  flow  case  with  axial  mean  flow  and  mean  shear,  Ux,2,  the  algebraic  equations  for  the 
fixed  points,  setting  D / Dtbij  =  0,  become 

(6n  +  l/3)(2  +  4bl2Sk/e)  -  4buSk/e  +  4 (Sk/e)  nil  /(4kS)  -  C,6„  -  2/3 
-\-C2{b\\b\\  +  612612  +  277/3)  —  8(fl/5)(5/s£)6i2  =  0 

612(2  +  4buSk/e)  -  2(622  +  l/3)Sfc/£  +  4(5fc/e)  II12  /(4^)  -  Ci612 
+^2612(611  +  622)  +  4(ft/ 5)(5fc/e)(6n  —  622)  =  0 

(622  +  l/3)(2  +  4buSk/e)  +  4(Sk/e)  ^22  /(4 kS)  -  Ctb22  -  2/3 
+C2(b22b22  +  612612  +  27  7 /3)  +  8(fi  /  5)(5  k  /  £)bx2  =  0 

(633  -I-  l/3)(2  +  4b\2Sk/e)  4-  4(Sk/e)Y\^3  f(4Sk)  —  C3633  —  2/3  +  £2(633633  +  27 7/3)  =  0 

The  last  equation  for  633  is  not  linearly  independent  as  6>j  =  0.  The  following  fixed  point  values  are 
taken  to  be  representative:  6f^  =  0.203,  6 ^  =  -0.156,  632  =  -0.143,  633  =  -0.06,  (Sk/e)^  =  5.54, 
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(PMco  =  1.73.  awe  inserted  into  the  algebraic  equations  above.  The  values  of  the  invariants 
corresponding  to  these  values  of  the  anisotropy  tensor  are:  //<„  =  -0.0604, 1 1 loo  =  0.0038,  Foo  = 
0.5613.  The  solution  of  the  algebraic  equations  describing  the  stationary  state  of  the  homogeneous 
shear  produces  the  following  values  for  the  calibration  coefficients, 

A%  =  -0.3388  -  0.06(i4g  +  Acw)  +  O.OlSCf0  -  0.008C£° 

Acn  =  - 14.35  -  2 Ac13  -  12.68^1  -  3.80/l|  -  5.0724fo  +  3.157C,00  +  0.5898Cf> 

A\ 2  =  -16.87  -  24^  +  28.7Ag  -  8.05A|  -  44.81A^0  -  12.17C°°  -  0.291C|°. 

The  model  has  been  left  as  general  as  possible  -  allowing  for  any  return  term  of  the  canonical  form, 
—C\ebij  +  C2 [*^  +  2///3£jj],  were  the  C,  are  not  necessarily  constants  but  have  achieved  their 
asymptotic  values  C,  =  Cf°.  Choosing  a  linear  Rotta  type  return  term  setting  the  nonlinear  return 
coefficient  to  zero,  C2  =  0,  produces 

/1|  =  -0.29  -  0.06(  Af0  -  Ac8) 

Acu  =  -3.6  +  5Af0  -  2Aci3  -  12.7A|  -  3.8A| 

Acu  =  -24.5  -  44.2Aj0  -  2^3  +  29ACS  -  SAC9 

where  the  following  values  of  the  free  parameters  have  been  chosen:  A3  =  0.8,  A9  =  —  1 .0,  A^0  = 
0.01,  A\3  —  0.  These  indicated  calibration  coefficients  are  appropriate  for  the  class  of  flows  in  which 
the  mean  shear  is  the  predominant  production  mechanism. 

The  component  form  of  the  rapid-pressure  used  in  the  equations  for  the  stationary  state  is 

Un  /VS*  =  C4612/3  +  C';\bnbn  +  6,2*12  +  2///3)6,2  +  Csbu(-W/S *) 

+[C6  +  C'l'(b\\  +  *22)1*11*12  +  C7(6,,  +  *22)6,2/3 

+Cs(*„  +  *22)*, 2(  W* / S*)  +  C9(^VS*)*,2(6,2*,2  -  *,,*22) 

IIS2  /2  ?2S*  =  C4*,2/3  +  <T(*22*22  +  *12*12  +  2/7/3)*,  2  +  C5bi2(W*/S*) 

+[C6  +  C4'(*ll  +  *22  )]*22*,2  +  CV(*ii  +  *22)6,2/3 

+C8(*ii  +  *22)*,2(W^*/S*)  +  C^-W*  /  S*)b\ 2(612*12  -  *,1*22) 

Uu  /2 92s*  =  1/2[C3  -  2//C3  +  3 77/Cf]  +  1/2 C4(*„  +  *22)  +  cn*,!  +  *22)6,2*12 
+£5(611  —  *22)(^*/2S*)  +  [C6  +  C^"(*„  +  *22)]*12*,2 
+  1/2C7(*11*11  -f  2*12*12  +  *22*22) 

+C8(*11*11  -  *22*22)(^*/2S')  +  C9(*11*22  “  *12*12)(*11  ~  *22)(W^*/2S*) 

1133  /2 q2S-  =  -2/3C4*,2  +  <T(*33*33  +  277/3)6,2  +  (C6  +  £$"(*„  +  *22))*12*33 
-2/3C7*l2(*ll  +  *22)- 


Appendix  3:  The  representation  of  the  integral  of  the  two-point  covariance 

The  general  form  of  the  fourth-order  tensor  polynomial  used  to  model  the  volume  integral  of  the 

derivative  of  the  two- point  covariance,  Xi}k j,  appearing  in  the  Reynolds  stress  equation  is 

Xtjki/  <  «p >=  M 6i}6ki  4-  A2(6ik6}i  4-  SnSjk) 

+  A36tJbki  +  A4bt:  Ski  +  As(bik6ji  4-  bn6jk  4-  6ikbji  +  6tib}k) 

+  Aebij^ki  +  Ajb^Ski  +  A&(b*k6ji  4-  b2t6}k  +  bikb2,  4-  b,tb2k) 

4-  A9bi}bki  4-  Aiofakbji  4-  btibjk) 

+  Anbitf,  4-  Aubfyki  4-  Al3{b}kb]l  4-  bfo*  4-  4- 

+  Aub2]b2kl  +  Als(b2kb2l+b2b2k) 

The  coefficients  with  all  free  parameters  set  to  zero,  as  derived  by  Ristorcelli  (1991)  or  in  Ristorcelli 
and  Lumley  (1987,1991),  are 

Ax  =  (111/7  4-  73)/27IId  —  F(420II  4-  239)/135//j 
A2  =  -(69//  4-  32)/27 Ild  +  F(420JI  +  257)/270//«( 

A3  =  (3 II  4-  4)/3 Ih  -  F(ll/10)/(1  4-  3//) 

A4  =  (15//  4-  1 1  )/3//<i  —  F(4/10)/(l  4-  3//) 

=  -3(1  4-  3II)/3IId  4-  F(3/10)/(l  4-  3//) 

A6=  -(102//4-61)/3//d 

/17  =  -2(33//  4-  20)/3 IId  -  F(6/10)/(l  4-  3//) 

A8  =  (42//  4-  23)/3IId 

A9  =  -(57 II  +  28)/3IId  -  F(3/10)/(l  4-  3//) 

A10  =  (15//  4-  U)/3IId  +  F(9/10)/(l  4-  31 1) 

An=  —(102//  4-  61)///j 
Ai2  =  -2(33//  4-  20)/IId 
A13  =  (42//  4-23)/ / /rf 

^H  =  0 

/t15=  0 

where  //rf  =  (1 4-3//)(74- 12//),  F  =  1 4-27/7/4-9//, where//  =  -1/2 and  ///  =  l/3bipbpjbji. 
The  rapid-pressure  correlation  integral,  Xpkj ,  appearing  in  the  heat-flux  equations,  is  used  to 
derive  the  >4,  through  the  joint- realizability  constraint.  The  general  form  of  the  third-order  tensor 
polynomial  used  to  model  Xpk j  is  included  for  completeness: 

Xvkj  =  D\  <  9up  >  6kj  4-  D2[<  9uk  >  6pj+  <  9uj  >  6P*] 

4-  D3<  9up  >  bkj  4-  /M<  9uk  >  bpj+  <  9u,  >  bpk] 

4-  £>5  <  9uv  >  b2kj  +  D%[<  9nk  >  6^4-  <  9uj  >  b2pk ] 

+  [Dibgpbkj  +  Ds(bqkbpj  4-  h?j^Pfc)]  <  9uq  > 

+  [Dvbqpbkj  4-  Dio{bqkbvj  4-  bq]  bpk )]  <  9uq  > 

+  [Dubqpb2}  +  D12(bqkb2pj  4-  bqj*pk)]  <  9uq  > 

+  [Dl3b2p6kj  4-  Du(b2k6PJ  +  b2j6pk)]  <  9uq  > 

+  [£>15^6*,  +  D^{b\kbp}  4-  b2q-bpk))  <  9uq  > 

+  ID* 7b2qpb2}  4-  Dia(b)kb2„  4-  b2qjb2pk)]  <  9uq  > 
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Where  the  D,  are  given  in  Ristorcelli  and  Lumley  (1987,1991)  as 

D\  =  — (312//2  +  149/7  -  21)/5//d  -  F(l/5)/(l  +  3//) 

D2  =  (48 II2  +  11-  14)/5 1  Id  +  F(3/10)/(l  +  3//) 

D3  =  -(324 II2  +  2221 1  +  17 )/IId  -  F( 3/10)/(l  +  3 //) 

£)4=  -3/(7+ 12//)  + F(9/10)/(l  +  3//) 

Ds=  -(102//  +  61)///^ 

£)6  =  (42//  +  23)///j 

£)7  =  -2(3//  +  4)/5//d  -  F(3/5)/(  1  +  3//) 

A*  =  27(2//  +  l)/5//ei 
£)9  =  (42//  +  23)///rf 
£>io  =  (42// +  23)///d 
£)13  =  -8(39//  +  22)/5 IId 
Du  =  2(24//  +  17)/5//d 
015=  -27/(1 +  3//) 

and  £>n  =  £>12  =  />i6  =  £>17  =  His  =  0.  The  first  term  of  each  of  the  coefficients  A{,  D,  satisfies 
the  constraints  of  continuity,  normalization,  realizability,  joint- realizability  and  2DMFI.  The  terms 
proportional  to  F  are  the  terms  that  contribute  to  the  rapid-pressure  when  the  turbulence  is  away 
from  both  the  realizable  and  the  2DMFI  limits. 

Appendix  4:  The  Cayley-Hamilton  theorem  generalization 

Reference  has  been  made  to  a  generalized  Cayley-Hamilton  theorem,  Rivlin  (1955),  relating  different 
powers  of  products  of  matrices. 

ABC  +ACB  +  BCA  +  BAC  +  CAB  +  CBA  = 

A{<  BC  >  -  <  B  ><  C  >)  +  B(<  CA  >  -  <  C  ><  A  >)  +  C(<  AB  >  -  <  A  ><  B  >) 
+(BC  +  CB)  <  A  >  +{CA  +  AC)  <  B  >  +(AB  +  BA)  <  C  > 

+  1[<  A  ><  B  ><  C  >  —  <  A  ><  BC  >  —  <  B  ><  AC  >  —  <  C  ><  AB  > 

+  <  ABC  >  +  <  CBA  >] 

The  Cayley-Hamilton  generalization  is  easily  derivable  from  the  Cayley-Hamilton  theorem  applied 
to  sums  and  differences  of  matrices  A,  B,C.  The  theorem  is  useful  in  eliminating  redundant  tensor 
bases  in  tensor  representation  theorems.  Here  <>  is  used  to  indicate  the  trace.  The  Cayley- 
Hamilton  theorem  for  the  anisotropy  tensor  is  b3  =  1/3  <  b3  >  1  +  1/2  <  62  >  6.  The  theorem  can 
also  be  used  to  express 

bSb  =  -[b2S  +  Sb2]+  <  bS  >  b  +  1/2  <  b2  >  S+  <  b2S  >  1 
bSb2  +  b2Sb  =  -1/3  <  b3  >  S+  <  b2S  >  b+  <  bS  >  b2 

Appendix  4:  The  rapid-pressure' models 

The  detailed  form  of  the  pressure-strain  models  referred  to  in  this  paper  are  as  follows: 

The  Launder,  Reece  &  Rodi  Model 

ntJ  =  -2C\ebij  +  ^ K SX]  +  C2K(bikSjk  +  bjkSik  -  ~h iSkiSij)  +  C3K(bikWjk  +  bjkW,k) 
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where 


The  Shih  &  Lumley  Model 


C,  =  1.5,  Ci  =  1.75,  C3  =  1.31 


n,;  =  -0£btJ  +  ^KStJ  +  12as K{bikS}k  +  b]kSlk~bklSk,6tJ)  +  1(2  -  7as)K(btkW}k  +  b]kWtk) 
+  -^K(b,tbtmSjm  +  bj[bimS,m  -  2bikSkibij-3bkiSkibij)  +  -K(b«bimW]rn  +  b}[bimWxm) 

where 

0  =  2+  ^  exp(-7.77/v//teT){72/  s/Ret  +  80.1  fn[l  +  62.4(-//  +  2.3///)]} 

F  =1+911  +  27///,  II  =  ~btJbt],  III  =  ±btjbjkbk, 


B  4  A'2  1  /  4  i\ 

Ret  _  _ —  <*5  =  tj:  (  1  +  7^*  )  • 

9  i/£  10  \  5  / 


The  Fu,  Launder  &  Tselepidakis  Model 

n,j  =  0i£b,j  +  02 e(bxkbkj  -  -bkibki6tj)+-KSij  +  1.2K(bxkSjk  +  bjkSik  -  -bkiSki6tJ) 

+  j^K(b,kW]k  +  bjkW,k)  +  -K(b,kbktS}t+b]kbkiS,i  -  2 bikSkibij  -  3bkiSkibi:>) 

+  \K(bikbklWjt  +  bjkbklW«)  -  jK[SII(b>kWjk  +  b]kWik)+\2(btkbklWlmbm]  +  bjkbklW,mbmi )] 


where 


0i  =  120//F1/2  +  2 F''2  -  2,  02  =  144 //F,/2 


The  Speziale,  Sarkar  &  Gatski  Model 


nu  =  -(2 Cxt  +  c;v)btJ  +  C2e  (bikbkj  -  l-bklbk,6i ,•)  +  (C3  -  c;ii>  )ks,} 
+C4K  (bikSjk  +  blkSik  -  | bklSk,6i j')  +  CsK{bikWjk  +  bjkWik ) 


where 


Ci  =  1.7,  Cj*  =  1.80,  C2  =  4.2 

C3  =  1  C3  =  1.30,  C4  =  1.25 
0 

C5  =  0.40,  I h  = 
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Figure  2c 


Figure  2d 
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Figure  5.  Time  evolution  for  the  anisotropy  tensor  in  plane  strain  for  e0/Sko  —  2.0.  The 

predictions  of  the  2DMFI  model -  ,  the  SL  model  — —  ,  and  the  direct  numerical 

simulations  of  Lee  and  Reynolds  (1985). 
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Figure  6.  Time  evolution  for  the  anisotropy  tensor  in  the  axisymmetric  contraction  for 

eo/Sko  =  0.179.  The  predictions  of  the  2DMFI  model -  ,  the  SL  model - ,  and  the 

direct  numerical  simulations  of  Lee  and  Reynolds  (1985). 
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